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SUMMARY 


A thermodynamic  model  of  combustion  and  radi- 
ation in  the  plume  of  a pyrotechnic  flare  has  been 
constructed  for  the  purpose  of  predicting  radiant 
efficiencies  of  untried  formulations. 

Plume  properties  and  flow  are  modeled  in  one 
dimension  only--paral 1 c 1 to  the  axis  of  the  plume. 
The  plume  composition  is  derived  from  the  equili- 
bration of  the  binder,  oxidizer,  admixed  air,  and 
a portion  of  the  metal  fuel  particles.  The  vapori- 
zation of  fuel  is  assumed  to  be  the  slow  step  in 
the  reaction.  Energy  losses  are  assumed  due  to 
radiaticn  only.  Energy  transfer  along  the  axis 
of  the  p« ume  by  radiative  or  other  mechanisms  is 
i gnored. 

The  starting  point  is  taken  in  the  gas  phase 
just  above  the  flare  surface.  All  of  the  binder 
and  oxidizer  are  assumed  to  have  reacted,  along 
with  enough  of  the  fuel  (5-10%)  to  raise  the  adia- 
batic flame  temperature  to  about  1000  K,  in  accord 
with  experimentally  deduced  vilues. 

Equilibrium  temperatures  and  species  concen- 
trations are  computed  at  closeiy.-spaced  points 
from  the  starting  point,  proceeding  parallel  to 
the  plume  axis,  until  well  into  the  tail  of  the 
plume  where  temperatures  have  dropped  to  the  point 
where  the  radiation  of  interest  is  essentially 
zero.  Increments  of  air  and  of  any  remaining  fuel 
are  added  to  each  point.  Radiation  losses  for  all 
of  the  important  radiative  processes  and  for  the 
processes  of  interest  are  computed  for  the  interval 
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between  each  pair  of  points  and  are,  subtracted 
from  the  enthalpy  used  in  calculating  equilibrium 
properties  at  the  next  point. 

For  a test  of  the  model  magnesium-fueled 
flares  v/ith  a sodium-containing  oxidizer  and  an 
organic  binder  were  chosen.  It  was  assumed  that 
the  most  important  radiative  processes  were  so- 
dium-D  emission  and  graybody  emission.  A radia- 
tive transfer  model  in  an  optically  dense  medium 
was  constructed  for  computing  sodium-D  emission. 
Graybody  radiation  was  taken  as  proportional  to 
the  fqurth  power  of  the  local  temperature. 

A computer  program  was  obtained  for  perform- 
ing the  calculations  by  modifying  a program  in 
the  literature. 

There  were  a number  of  parameters  whose  values 
could  not  be  determined  a priori.  These  included 
the  rates  of  air  mixing  and  of  fuel  consumption, 
the  absolute  rates  of  radiation,  and  the  plume 
spreading  angle.  The  functional  form  of  these 
quantities  was  deduced  from  theory  or  arbitrarily, 
and  the  absolute  values  were  adjusted  so  as  to 
give  agreement  with  observations  of  plume  shape, 
the  point  of  maximum  brightness,  the  total  visible 
length,  and  the  luminous  efficiency  for  a limited 
set  of  observations. 

Computations  were  then  performed  for  a wide 
range  of  fuel/oxidizer  ratios,  oxidizer  formulas, 
binder  percentages,  ambient  pressures,  and  am- 
bient compositions  and  compared  with  experimen- 
tal values  from  the  literature.  Agreement 
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between  computed  and  literature  values  of  lu- 
minous efficiency  varied  from  moderate  to  good. 

Modifications  were  indicated  by  which  the  pro- 
gram can  be  used  for  formulas  where  other  types  of 
radiation  are  of  interest.  An  example  suggested  Q 
is  the  calculation  of  the  green  radiation  at  5290  A 
from  the  En  excited  electronic  state  of  the  BaCl 
radi cal . 

It  was  recommended  that  the  necessary  data  be 
acquired  and  .modi fi cations  made  so  that  the  program 
could  be  applied  to  colored  flares  or  to  infrared 
flares.  It  was  also  recommended  that  the  model 
be  refined  (1)  by  a reoptimization  of  the  arbitrary 
input  parameters  (rate  of  air  mixing,  etc.),  (2)  by 
modifying  the  fuel  input  (vaporization)  to  take  ac- 
count of  local  oxygen  concentration,  and  (3)  by 
including  a consideration  of  the  heat  feedback 
parallel  to  the  plume  axis. 


PREFACE 


This  work  is  part  of  a' larger  program  directed 
at  a better  understanding  of  pyrotechnic  flcire  com- 
bustion and  radiation.  Besides  modeling,  the  pro- 
gram includes  high  resolution  spectroscopy  and  band 
radiometry  in  the  visible  and  infrared  wavelengths. 

Thanks  are  due  to  Drs.  Harold  Sabbagh,  Henry 
Webster,  Carl  Dinerman,  and  John  O'Benar  of  NWSC 
for  helpful  comments  on  the  manuscript. 

The  typing  and  editing  of  this  report  by 
Miss  Karla  Harker  is  greatly  appreciated. 
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INTRODUCTION 


Many  experimental  and  theoretical  studies  of 
the  combustion  of  flares  have  been  made  to  determine 
the  details  of  the  combustion  and  light  emission 
processes  so  that  a better  understanding  of  their 
operation  could  point  the  way  to  devising  flares 
of  greater  luminosity  or  color  purity.1-7  Model- 
ing of  related  phenomena  has  also  been  done.8 

The  task  of  the  present  study  has  been  to  con- 
struct a theoretical  model  of  flare  combustion 
which  would  incorporate  all  the  experimental  ob- 
servations and  at  the  same  time  test  them  for  con- 
sistency. The  goal  is  to  be  able  to  predict 
radiant  efficiency  of  new  formulations  at  selected 
wavelengths  or  spectral  bands.  In  addition,  it 
may  be  possible  to  gain  additional  information  on 
details  of  the  flare  combusti on--s uch  as  actual 
flame  temperatures,  the  emissivity  of  the  smoke, 
and  burning  times  of  the  fuel  particles. 

Some  of  this  work  has  been  reported  previ- 
ously. 0 
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EXPERIMENTAL  BACKGROUND 


For  initial  tests  of  the  model,  we  have  taken 
the  magnesium/sodium  nitrate/organic  binder  flare 
composition  as  an  example,  since  this  composition 
represents  the  most  common  military  illuminating 
flare,  and  the  one  most  thoroughly  observed. 

Various  evidence  indicates  a temperature  of  about 
1000  K at  the  burning  surface:  It  has  been  ob- 

served,7 using  high  speed  photography,  that  most 
magnesium  particles  (m.p.  922  K)  melt  and  ignite 
a few  milliseconds  before  ejection  from  the  solid, 
while  some  ignite  a short  distance  away  from  the 
surface.  Furthermore,  decomposition  of  sodium 
nitrate,  which  produces  the  gas  to  blow  the  mag- 
nesium out  of  the  flare,  has  been  shown  to  take 
place  very  rapidly  at  1000  K10  at  perhaps  the  same 
rate  as  in  the  surface  of  a burning  flare. 

Finally,  ignition  temperatures  of  magnesium  dust 
in  air  have  been  reported  to  be  about  900  K.11 

Plume  temperatures  have  been  determined  by 
thermocouples12  and  by  absolute  brightness  measure- 
ments13 at  the  frequencies  of  the  two  maxima  in 

O 

the  broadened  sodium  D line  K5800  and  6000  A). 

Both  types  of  measurements  yield  maximum  tempera- 
tures in  the  vicinity  of  the  commuted  adiabatic 
temperature,  although  they  are  uncertain  by  several 
hundred  degrees.  The  position  of  the  maximum 
temperature  is  found  within  a few  inches  of  the 
surface  in  the  thermocouple  measurements  cited. 
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Observations14  indicate  that  the  region  of 
greatest  brightness  begins  about  one  diameter  from 
the  surface.  For  stoichiometric  mixes,  the  brightest 
region  is  fairly  concentrated,  whereas  for  fuel- 
ri ch  mixes  it  is  well  spread  out  along  the  length 
of  the  plume.3 

The  shape  of  the  flame  and  smoke  plume  de- 
pends on  the  orientation  of  the  flare  and  on  the 
surrounding  air  flow.  For  a flare  burning  upright 
in  still  air,  the  plume  can  be  approximated  as 
cone-shaped,  with  an  apex  angle  of  15-30°.  Plumes 
of  large  diameter  (11  cm)  flares  may  be  luminous 
for  a length  of  1.5  to  2 meters  to  a point  where 
the  photometric  brightness  has  decreased  by  a fac- 
tor of  about  106  from  the  maximum  value. 

Power  spectra  indicate  that  most  of  the  lumi- 
nous intensity  is  from  the  sodium  D radiation,1'2 
a doublet  of  wavelength  589  nm.  This  radiation 
accounts  for  about  8%  of  the  total  heat  of  combus- 
tion in  flares  which  have  been  optimized  for 
efficiency  of  luminous  output. 
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THE  MODEL 
Flow  in  the  Plume 


Combustion  of  a flare  is  quite  complicated.  In 
the  body  of  the  flare  many  reactions  occur  as  the 
binder  and  oxidizer  break  down.  Three  phases  are 
present  near  the  surface  as  bubbles  of  gas  push 
their  way  past  remaining  liquid  and  solid. 

In  the  plume  itself,  the  overall  flow  of  material 
in  three  dimensions  is  complicated  by  strong  local 
turbulence.  A further  complication  of  the  flow  is 
a detectably  different  movement  of  the  gases  and  of 
the  heavier  particles.15 

Equilibration  of  the  reactants  is  slowed  in 
many  formulations  by  the  high  reaction  temperatures 
of  the  metal  particles. 

Transfer  of  heat  within  and  out  of  the  plume 
occurs  by  a variety  of  radiation  and  conduction 
processes.  The  radiation  is  complicated  in  many 
flares  by  the  light  scattering  of  a highly  reflec- 
tive smoke. 

The  primary  goal  of  this  project  was  to  pre- 
dict radiative  output  rather  than  the  details  of 
combustion  and  flow.  Therefore  many  simplifying 
assumptions  were  made  regarding  these  latter 
phenomena . 

First  of  all,  the  plume  was  modeled  in  one 
dimension  with  the  flow  of  gases  in  one  direc- 
tion only--parallel  to  the  plume  axis.  This  most 
closely  approximates  the  case  of  a flare  burning 
upright,  or  else  pointed  downward  in  a forced 
downdraft . 
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The  expansion  of  the  plume  perpendicular  to 
the  axis  depends  not  only  on  the  amount  of  air 
entrained  but  on  viscous  friction  at  the  plume 
surface.  Since  this  friction  was  not  included  in 
the  model,  the  plume  shape  was  arbitrarily  speci- 
fied as  a cone,  with  apex  behind  the  burning  sur- 
face, as  shown  in  Figure  1.  Since  the  properties 
perpendicular  to  the  axis  were  taken  as  constant 
the  problem  is  still  one-dimensional. 

While  the  turbulence  of  the  plume  flow  was 
not  modeled,  one  of  the  major  effects  of  turbu- 
lence--the  entrainment  of  the  surrounding  environ- 
ment (usually  air)--was  included.  The  model  as- 
sumes continuous  addition  of  air  along  the  plume 
at  an  arbitrarily  specified  rate. 

Conditions  for  the  conservation  of  momentum 
in  the  plume  are  not  imposed  by  the  model.  How- 
ever an  ad  hoc  calculation  is  made  of  the  momen- 
tum change  in  the  forward  direction  due  to  addi- 
tion of  environment  (air). 

We  start  with  mv  = mF/Ap^,  where  m is  the  mass 
of  the  cross-sectional  slice  of  plume  under  con- 
sideration, F is  the  total  mass  flow  across  a 
cross-sectional  plane,  A is  the  cross-sectional 
area,  and  p,  is  the  plume  density.  The  change  in 
plume  density  due  only  to  the  addition  of  the 
environment  is  given  approximately  by 


P1(P3  I P1} 
dir  - P z[i  + r)  ■ » 
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1.  Plume  parameters,  labeled  with  the  variable  names 
the  computer  program.  See  Table  A-l  in  Appendix  A. 


where  pg  is  the  air  density  and  R is  the  weight 
ratio  of  air  to  initial  composition. 

The  change  in  A is  computed  from  the  assumed 
plume  shape;  and  the  change  in  the  product  mF. equals 
the  initial  value,  mQFo,  times  R.  We  then  find: 


mixing 


m F (I  + R) 

ps  + p,  dR 

A p 

P2  dx 

d In  A 
<J  + r)-3x- 


(1) 


where  x Is  parallel  to  the  plume  axis. 

If  at  any  time  the  quantity  in  brackets  is  found 
to  be  positive  a warning  is  printed  out.  When 
a cone-shaped  plume  is  assumed  it  can  be  shown 
that  the  admission  of  ai r at  a rate  constant  in 
distance  will  generally  cause  this  quantity  to 
be  less  than  zero,  as  desired. 


Temperature  and  Energy 

A complete  model  of  a flare  flame  in  one  di- 
mension would  satisfy  the  condition  that  the  net 
energy  flow  into  a point  is  zero  (compare  Eqn.  5.2 
of  Ref.  16) : 


37(ki7>  ' cpV&  + «<T>  - R<T>  ’ 0 U) 

The  first  term  represents  heat  transport  in  the 
direction  of  the  flow  by  simple  conductivity,  radi- 
ation, and  local  turbulence.  The  second  term  repre- 
sents net  heat  transport  due  to  flow  of  material 
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containing  a temperature  gradient.  The  third 
term  represents  generation  of  heat  by  chemical 
reaction.  The  fourth  term  is  loss  of  heat  to  the 
outside  (perpendicular  to  the  direction  of  flow) 
by  radiation  and  simple  conduction.  Q and  R are 
strongly  dependent  on  temperature.  In  general,  the 
thermal  conductivity  k,  heat  capacity  c^,  velocity 
v , and  density  p have  a much  weaker  temperature 

A 

dependence . 

With  relatively  minor  simplifications,  Eqn.  (2) 
could  reasonably  be  applied  to  the  case  of  a flare 
plume,  using  the  numerical  techniques  employed  here 
However,  for  a first  effort,  to  save  considerable 
programming  and  computer  time,  the  first  term  in 
this  equation  (the  thermal  conductivity)  was 
dropped.  Since  this  term  represents  the  thermal 
feedback,  whi ch  aids  greatly  in  driving  the  reaction 
particularly  in  the  early  stages,  it  was  neces- 
sary to  arbitrarily  specify  Q rather  than  model 
it  on  physical  principles,  in  order  to  get  realis- 
tic results.  This  was  done  by  arbitrarily  speci- 
fying the  rate  at  which  the  metal  fuel  enters  the 
chemical  reaction  by  vaporization,  as  will  be  ex- 
plained in  the  section  "Reactions  in  the  Plume". 

The  loss  of  heat  to  the  outside,  represented 
by  R(T),  was  broken  into  several  parts  to  repre- 
sent different  types  of  radiation  processes.  The 
functional  form  of  each  was  modeled  to  a first 
approximation.  Non-radi ati ve  losses  to  the  out- 
side can  be  safely  ignored.  Thermal  conductivity 
by  molecular  diffusion  should  b,e  relatively  small. 
Heat  loss  due  to  mixing  in  of  air  at  the  plume 
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surface  is  already  taken  into  account  by  including 
this  air  with  the  reactants. 

Equation  (2)  without  the  first  term  is  solved 
numerically  using  a modification  of  the  NASA  com- 
puter program  for  complex  equilibria.17  The  origi- 
nal program  computes  equilibrium  compositions  and 
temperatures  from  a given  set  of  reactant  materials 
and  their  enthalpies.  Modifications  were  added 
to  cause  the  program  to  automatically  recompute 
equilibrium  a specified  number  of  times  to  repre- 
sent closely  spaced  points  along  the  plume.  At  each 
point,  increments  of  air  and  of  magnesium  are 
added  to  the  amount  of  initial  reactants,  to  repre- 
sent the  continual  mixing  in  of  air  and  vaporiza- 
tion of  magnesium.  Radiation  losses  are  computed 
based  on  temperature  and  emitter  (sodium  atom) 
concentrations  at  each  point.  These  losses  are 
subtracted  from  the  enthalpy  used  in  computing 
equilibrium  at  the  next  point.  In  this  way,  a 
temperature,  composition,  and  radiation  profile 
are  constructed  along  the  length  of  the  plume.  The 
radiant  outputs  are  added  up  for  all  of  the  points, 
separately  for  each  emission  process,  to  obtain  the 
total  output  of  the  amount  of  material  considered. 

The  values  of  p computed  thermodynamically  at 
each  point  are  used  to  recompute  v , and  thus  vari- 

A 

ations  of  c , v , and  p with  changes  in  composition, 
P x 

temperature,  and  degree  of  reaction  are  automati- 
cally taken  into  account. 

If  it  were  desired  to  include  thermal  feedback, 
term  number  1 of  tqn.  (2),  several  iterations  would 
need  to  be  made  at  each  distance  increment  in  the 
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plume.  The  value  of  k which  includes  radiation, 
conduction,  and  turbulence  contributions,  would 
have  to  be  specified,  along  with  its  temperature 
dependence,  perhaps  as  one  or  two  adjustable  para- 
meters . 


Initial  Conditions 

The  ideal  starting  point  for  the  computations 
would  be  in  the  flare  body  at  ambient  temperature, 
well  behind  the  burning  surface.  However,  the  il- 
lumination produced  by  the  flare  body  is  negli- 
gible, and  probably  its  radiation  at  longer  wave- 
lengths can  also  be  safely  ignored  because  of  its 
small  size  and  low  temperature  compared  to  the 
plume.  Therefore,  to  avoid  the  complications  of 
modeling  the  processes  taking  place  in  the  solid, 
particularly  at  the  surface,  and  having  to  con- 
sider heat  transfer  to  the  surface,  the  computa- 
tion is  started  in  the  gas  phase  at  the  burning 
surface.  For  the  magnesium/sodium  nitrate  flare, 
a temperature  of  1000  K is  assumed.  This  is 
in  reasonable  agreement  with  the  observations 
quoted  earlier. 

We  must  require  that  the  change  from  ambient 
temperature  to  1000  K occurs  adi abati cal ly . This 
is  because  of  the  neglect  of  the  heat  feedback 
(term  number  1 of  Eqn.  (2))  and  the  assumption  of 
no  heat  loss  by  the  flare  body.  Just  enough  of  the 
original  fuel  is  assumed  to  have  vaporized  and  en- 
tered the  reaction  to  fulfill  this  adiabatic 
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condition.  The  remainder  of  the  fuel  has  not  yet 
entered  the  reaction.  Under  equilibrium  conditions 
at  1000  K most  oxidizers  and  binders  are  well  de- 
composed. 


Reactions  in  the  Plume 

t 

Binder  and  oxidizer  ars  assumed  to  have  com- 
pletely decomposed  on  leaving  the  flare  surface. 

The  magnesium  is  assumed  to  leave  the  surface 
in  the  form  of  the  original  particles,  which  grad- 
ually vaporize  during  flight. 

Mixing  and  chemical  equilibration  among  the 
decomposition  products  of  binder  and  oxidizer, 
the  magnesium  vapor,  and  the  entrained  air  is  as- 
sumed to  be  instantaneous.  This  is  in  agreement 
with  propellant  models,  where  the  thickness  of  the 
flame  reaction  zone  above  the  burning  surface  has 
been  calculated  to  be  of  the  order  of  millimeters.8 

The  rate  of  reaction  is  then  determined  by  the 
rate  of  vaporization  of  the  magnesium  and  the  rate 
of  entrainment  of  air.  An  expression  for  the  rate 
of  vaporization  of  the  magnesium  particles  as  a 
function  of  local  temperature  and  chemical  en- 
vironment (oxygen  concentration)  based  on  avail- 
able particle  burning  models  could  be  incorpo- 
rated into  the  flare  combustion  model.  But  be- 
cause of  the  mutual  interaction  between  temperature 
and  rate  of  burning,  it  would  then  be  necessary  to 
include  the  heat  feedback  term  of  Eqn.  (2)  in  order 
to  get  the  reaction  started  in  a reasonable  time. 
Since  this  term  was  omitted,  the  next  best  alterna- 
tive was  to  specify  a rate  of  vaporization  of 
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magnesium.  For  want  of  better  knowledge,  we  as- 
sumed that  the  magnesium  vaporized  at  a constant 
rate  In  time  until  consumed.  In  fact,  this  Is  a 
crude  compromise  between  the  acceleration  of  the 
burning  rate  due  to  the  increase  in  temperature 
when  proceeding  out  into  the  plume  and  the  decele- 
rating effects  of  the  decrease  in  particle  size 
(surface  area)  and  the  accumulation  of  combustion 
products  near  the  surface. 

It  should  be  pointed  out,  however,  that  due  to 
neglect  of  heat  feedback,  the  computed  progress 
of  the  temperature  of  reaction  will  lag  behind 
the  real  temperature  for  any  degree  of  combustion. 

An  option  was  included  in  the  computation  to 
terminate  vaporization  of  the  magnesium  before  the 
particles  are  completely  consumed.  This  option  is 
useful  for  testing  the  hypothesis  that  certain 
fuels  (e.g.  aluminum)  give  lower  output  than  ex- 
pected due  to  incomplete  combustion. 

The  magnesium  particles  were  assumed  to  be  at 
the  local  temperature  or  at  the  normal  boiling 
point,  whichever  was  lower.  Appropriate  correc- 
tions to  the  total  enthalpy  were  made  for  the 
heat  content  of  the  unvaporized  magnesium. 

Radi ati on 

As  stated  before,  energy  loss  to  the  outside 
is  taken  to  consis.  entirely  of  radiation.  A, num- 
ber of  kinds  of  radiation  processes  rnay  be  included, 
each  with  its  own  dependence  on  temperature  and 
concentration  of  emitting  species.  Possibilities 
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include  simple  atomic  line  emission,  atomic  line 
emission  in  an  optically  dense  medium,  molecular 
band  emission,  and  graybody  radiation. 

The  computer  program  is  presently  set  up  with 
two  radiation  terms--sodium-D  line  emission  in  an 
optically  dense  medium  and  graybody  radiation. 
Thus,  the  program  can  compute  luminous  intensity 
for  any  case  where  emission  of  a molecular  species 
or  of  an  atomic  species  other  than  sodium  does  not 
dissipate  a significant  {>2%)  portion  of  the  en- 
thalpy. For  somewhat  greater  accuracy,  the  emis- 
sion of  the  sooiu.Ti  doublets  at  568  nm  and  at  616  nm 
could  have  been  included. 

Sodium  D Line  Emission 

In  the  usual  case  where  the  concentration  of 
emitters  is  low  enough  such  that  reabsorption  is 
negligible,  the  increase  in  intensity  of  a beam 
traveling  from  one  side  of  the  plume  to  the  other 
is  given  by 

3"r  = ?X7  exP(~hc/xkT) , 


is  Planck's  constant 
is  Boltzmann's  constant 
is  the  velocity  of  light 
is  the  wavelength 


where 

h 

k 

c 

A 
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t is  the  natural  emission  lifetime,  and 
N is  the  number  of  emitting  atoms  or 
molecules  per  unit  volume. 

The  total  emerging  intensity  is  the  integral  of 
this  expression  across  the  plume;  i.e.  over  r. 

When  N/t  is  large,  or  the  path  is  long,  re- 
absorption must  also  be  considered.  Since  absorp- 
tion is  dependent  on  the  beam  intensity  and  on  the 
probability  of  absorption,  both  of  which  are  fre- 
quency dependent,  the  intensity  must  be  computed 
over  small  frequency  intervals.  Including  absorp- 
tion, we  have  (compare  Eqn.  25,  Ref.  18); 


dI  A _ 

w - 


B 


A 1 dN 
c d.\ 


TTT  iT  exp(-hc/xkT) 


(3) 


where  I represents  the  intensity  in  unit  wave- 
length interval  centered  at  A,  B = A3/2hcT  is  the 
Einstein  transition  probability  per  unit  frequency 
for  absorption,  and  dN/dA  is  the  number  density 
of  emitting  or  absorbing  atoms  per  unit  wavelength. 
We  rewrite  this  as 


dr 


= “kA^A 


IA(eq)], 


(4) 
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where 


= A4  dN 
x FttcT  <3T 


is  the  absorption  coefficient  at  wavelength  X,  and 
I x ( eq ) is  the  beam  intensity  at  thermodynamic 
equilibrium  in  the  wavelength  interval  x to  x + dx. 

The  total  intensity  1^  emerging  from  the  flare 
plume  will  be  approximately  the  integral  of  Eqn.  (4) 
across  the  plume  (alonq  r)  provided  that  vari- 
ations in  temperature  and  emitter  density  along 
the  plume  length  parallel  to  x are  not  too  abrupt. 

In  a real  flare  the  light  path  traverses  re- 
gions with  temperatures  ranging  from  very  high  to 
ambient,  which  causes  an  absorption  dip  at  the  cen- 
ter of  the  spectrum  of  the  line.  For  the  purpose 
of  including  this  feature  in  our  calculation  of 
the  emerging  Na  D radiation  we  will  make  an  ex- 
ception to  our  basic  assumption  of  uniform  proper- 
ties across  the  plume.  We  approximate  the  real 
situation  by  assuming  that  surrounding  the  plume, 
whose  local  diameter  is  1 , ir  a thin  sheath  of 
thickness  lg,  and  zero  temperature,  and  which  has 
the  same  sodium  atom  concentration  as  within  the 
plume.  (Light  emission  at  ambient  temperature  is 
negligible.  The  use  of  zero  temperature  is  alge- 
braically much  simpler).  Light  reabsorbed  here  will 
be  considered  as  not  lost  to  the  plume  enthalpy. 
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Equation  (4)  then  integrates  to 


I 


x 


Ix(eq)  axp(-klslg)(l  - expi-k^l,)) 


(5) 


Subscripts  1 and  2 refer  to  the  hot  and  cold  paths, 
respecti vely . 

To  get  the  total  intensity  we  must  integrate 
Eqn.  (5)  over  all  wavelengths.  To  do  this  we  need 
dN/dx,  in  the  expression  for  k^,  i.e.  we  need  to 
know  the  line  shape. 

Although  the  line  shape  is  strongly  influenced 
by  both  collision  broadening  and  Doppler  broaden- 
ing, the  effect  of  the  latter  becomes  negligible 
in  the  wings  of  the  line,  which  is  where  nearly 
all  the  radiation  emerges  for  the  high  sodium 
atom  densities  typical  of  illuminating  flares. 

We  may  then  assume  a Lorentzian  line  shape,  so  that 


dN  s N 

tfX  TT 


2/  A X 

1 + L2(x  -xo)/axlJs 


where  ax  is  the  half-width  at  half  height  for  the 

L 

line  as  broadened  by  collisions,  and  x is  the  wave- 

o 

length  at  the  line  center. 

On  substituting  this  formula  for  dN/ dx  into  kA 
and  in  turn  into  Eqn.  (5),  an  expression  was  ob- 
tained which  could  not  be  integrated  in  closed  or 
in  series  form.  A numerical  integration  would  have 
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been  possible,  but  since  A A , 1 and  1 vary  at 
each  point  along  the  combustion  path,  the  integra- 
tion would  have  to  be  repeated  many  times,  and  would 
increase  the  costs  of  computation. 

Since  only  the  far  wings  of  the  line  contri- 
bute significantly  to  the  emerging  light  intensity 
we  can  make  some  approximations: 

(1)  Combining  the  above  expressions  for  kx 
and  dN/dA  we  obtain 


kx  = k0/U  + (2(x  - Ao)/AAl)]2. 


kQ  = A4N/4tt2ctAAl  is  nearly  constant  over  the  line 
and  equals  che  absorption  coefficient  at  the  center 
of  the  1 ine  when  A = Aq. 

For  mathematical  simplicity  we  use  the  same 
values  of  aa  and  N (and  hence  of  k ) in  the  hot 

L A 

and  cold  regions.  This  approximation  is  not 
serious  in  regard  to  AAl  since  it  is  only  weakly 
dependent  on  temperature,  because  the  line  is  pres- 
sure broadened.  N varies  much  more  drastically. 

On  the  other  hand,  the  cold  portion  of  the  path 
must  be  relatively  longer  lower  temperatures. 
Therefore  we  will  consider  1 and  1 as  being  ef- 
fective  lengths  rather'than  actual  distances. 

Since  we  will  use  the  values  of  AAu  and  N 
computed  from  the  high  temperature  properties,  the 
error  will  be  in  the  estimation  of  the  absorption 
dip.  Since  this  dip  is  only  about  20%  of  the  total 


intensity,  errors  in  its  calculation  do  not  seri- 
ously affect  the  estimation  of  total  output. 

(2)  It  can  be  shown  from  Eqn.  (5)  that  there 
is  a maximum  in  I A as  a function  of  and  that 
the  wavelength  xm  at  this  point  is  given  by 


AX  / k 1 \% 

Am  ' Ao  = 1 Ty  n (T  /"l  ~ '+“77'  " 


In  evaluating  the  right  side  of  this  equation  we 
can,  to  a good  first  approximation,  set  x = Xq 
in  the  expression  given  earlier  for  kQ. 

In  the  limit 


k 1 » Y- -vX°  1 n ( 1 /I  + l) 

O 1 X - X 12 


we  have 


-)■ 


Tne  inequality  expresses  the  condition  that  the 
optical  thickness  at  line  center  should  be  much 
greater  than  10,  where  unit  optical  thickness  is 
defined  as  attenuation  by  the  factor  e.  This  should 
be  true  for  most  cases  of  practical  interest. 
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We  let  r = ] /l  and  y = (x  - x )/(x  - x ). 

12*  O III  O 

Substituting  into  the  expression  for  k^  and  into 
Eqn.  (5)  and  integrating  over  X we  obtain 


1 - Ix(eq)t», 


X0>/“  (r  +ir/r*' 

•/-  00 

x[2  - (r  + 3)"l/yS]dy 


(6) 


(3)  Since  the  integral  depends  only  on  r, 
if  we  assume  a constant  ratio  of  hot  path  length 
to  cold  path  length  in  the  plume,  we  need  perform 
the  integration  only  once. 

Since  we  assume  that  the  line  shape  is  domina 
ted  by  collision  broadening  we  evaluate  the  line 
width  by: 


AX,  = ( 2 X2o  2N_  /ttc)  [2trRT ( 2/M  + 2/M  ) f 

(compare  Eqn.  (104)  of  Ref.  18),  where 

aL 2 is  the  cross  section  for  transition  due  to 
col  1 i si ons  , 

N is  the  number  density  of  all  molecules, 

Mt  is  the  atomic  weight  of  the  emitting  atom 
or  molecule 

M is  the  average  molecular  weight  of  the  gases 
(in  the  plume) 

Taking  oL2  = 50  x 10"16  cm2  (Ref.  18  p.  171), 

M.  = 30,  t = 10"8  sec,  X = 5.89  x 10'5  cm,  N = 
Avogadro's  number  x PV/RT,  and  using  the  known 


19 


values  of  h,  c,  k,  R,  and  we  obtain 

I,(eq)(xm  - X ) = 32 2 [exp (-2441 6/T)][NPl  //Tln(r  + 2)]* 

Amo  i 

in  ergs  sec-1  cm"3*  Thus  the  total  intensity  may 
be  evaluated  from  Eqn.  (6). 

It  is  interesting  to  note  that  with  the  approxi- 
mations used, the  total  emerging  intensity  is  pro- 
portional to  the  square  root  of  sodium  atom  concen- 
tration (N)  and  path  length  (1  ),  rather  than  to 
the  first  power  of  these  variables  as  in  the  op- 
tically thin  case. 

The  simplifications  made  in  deriving  the 
radiative  transfer  equations  of  this  section  are 
based  on  the  assumption  of  a high  optical  den- 
sity--as  stated  earlier,  an  optical  density  much 
greater  than  10  at  line  center.  We  may  con- 
veniently estimate  these  optical  densities  from  the 
depth  of  the  absorption  dip  in  the  spectrum  of  the 
broadened  line.  An  examination  of  broadened 
sodium  D spectra  taken  by  Douda19  over  a variety 
of  ambient  pressures  and  sodium  concentrations 
indicates  that  this  treatment  will  not  hold  well 
enough  for  ambient  pressures  less  than  about  30 
Torr  at  high  sodium  content  or  150  Torr  at  a 5 
percent  NaN0&  content. 

Graybody  Emission 

As  is  well  known,  graybody  radiation  is 
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proportional  to  the  area  of  emitting  surface,  the 
surface  emissivity,  and  the  fourth  power  of  the 
temperature . 

In  a smoke,  surface  area  may  refer  to  either 
the  total  area  of  all  the  particles  or  the  area 
of  the  periphery  of  the  cloud,  depending  on  whe- 
ther the  cloud  is  optically  thin  or  thick,  respec- 
tively. Computations  indicate  that  typical 
magnesium  oxide  smoke  clouds  are  highly  reflective 
and  therefore  may  be  considered  as  optically  thin. 

Since  we  lacked  sufficient  information  about 
particle  sizes  or  the  emissivity  (and  its  tem- 
perature dependence)  of  the  magnesium  oxide  pro- 
duced we  have  expressed  the  time  rate  of  gray- 
body  radiation  loss  per  unit  weight  of  original 
composition20  by  the  simple  product  aj4,  where  a 
is  a constant  to  be  determined  by  trial,  and'is 
valid  whenever  magnesium  oxide  is  the  principal 
condensed  phase  produced. 

A further  calculation  compares  this  radiation 
with  the  blackbody  radiation  which  a flare  plume 
of  this  size  would  emit  (if  the  smoke  were  op- 
tically dense).  The  ratio  is  generally  about 
0.1.  It  can  be  thought  of  as  the  "emissivity" 
of  the  plume  boundary. 

The  luminous  intensity  of  the  graybody  radi- 
ation is  computed  and  presented  separately  from 
the  luminous  intensity  of  the  sodium  D radiation. 

An  option  is  also  included  for  computing  the 
portion  of  the’ graybody  radiation  which  is"  emitted 
into  selected  wavelength  bands  of  the  spectrum. 
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COMPARISON  WITH  EXPERIMENT 


Before  numerical  calculations  can  be  made, 
values  must  be  assigned  to  various  parameters 
which  are  included  in  the  model,  but  which  are  not 
known  a priori  and  were  not  derived  theoretically. 
As  presently  set  up,  this  includes  the  plume 
shape,  the  burning  rate  (surface  recession),  the 
proportion  of  admixed  air,  the  vaporization  rate 
of  metal  fuel,  adjustment  factors  to  the  semi- 
theoretical  expressions  for  the  radiation  proces- 
ses, and  the  initial  temperature. 

In  the  following  pages  it  is  described  how 
this  was  done  using  data  from  flares  of  the 
Mg/NaNO  /binder  composition.  Comparisons  are 
then  presented  of  calculated  and  experimental 
results  of  this  type  of  formula  with  a range  of 
fuel/oxidizer  ratios,  environmental  pressures, 
and  environmental  compositions,  and  with  a small 
range  of  binder  content.  A comparison  of  pre- 
dicted and  experimental  results  for  a selection 
of  oxidizer  compounds  is  also  presented. 

Adjustment  of  Parameters 

The  values  of  the  adjustable  parameters  de- 
pend in  a complex  manner  on  the  experimental 
observations  which  are  to  be  satisfied.  For 
instance,  the  graybody  constant,  the  plume  angle, 
and  the  rate  of  air  mixing  together  determine  the 
length  of  the  visible  part  of  the  plume.  The 
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particle  vaporization  rate,  the  plume  angle,  and 
the  rate  of  air  mixing  are  all  important  in  de- 
termining the  position  of  the  brightest  spot.  The 
graybody  and  sodium  D radiation  constants,  as  well 
as  the  proportion  of  admixed  air  strongly  affect 
the  predicted  luminous  output. 

Because  of  the  strongly  nonlinear  effect  of 
the  model  parameters  on  the  computed  results,  no 
systematic  procedure  was  found  for  adjusting  the 
parameters  to  optimize  the  fit  of  predicted  values 
to  observed  values.  However,  a large  number  of 
computations  were  performed  using  the  formula  58% 
magnesium,  37.5%  sodium  nitrate,  4.5%  epoxy  binder 
and  a 10.8  cm  diameter,  and  a best  set  of  para- 
meters was  selected. 

The  plume  shape  was  assumed  for  simplicity 
to  be  conical.  This  is  also  in  reasonable  agree- 
ment with  observations  of  flares  burning  upright. 

A plume  angle  of  0.6  radians  was  chosen.  This  is 
at  the  upper  limit  of  observed  values.  The  use 
of  a narrower  angle  causes  too  long  a visible 
plume  with  reasonable  values  of  the  sodium  ra- 
diation constant  and  of  the  proportion  of  admixed 
air. 

A burning  rate  of  0.25  cm/sec  similar  to  that 
of  the  Mk  45  Aircraft  Parachute  Flare  was  assumed. 
The  computed  luminous  output  as  a function  of 
burning  rate  (only)  passes  through  a broad  maxi- 
mum at  a value  of  about  0.5  cm/sec,  so  that  the 
results  are  moderately  sensitive  to  decreases  in 
burning  rate  but  rather  insensitive  to  assumed 
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The  initial  temperature  above  the  burning  sur- 
face is  assumed  to  be  1000  K.  A small  portion  of 
the  magnesium  is  assumed  to  have  entered  the 
reaction  at  this  point,  so  that  a computation  of 
the  initial,  adiabatic  temperature  yields  this  num- 
ber. However,  the  results  of  luminous  output  are 
affected  very  little  even  by  a value  several  hun- 
dred degrees  either  way.  This  is  because  so  little 
luminous  output  occurs  at  the  very  short  time  the 
temperature  remains  at  this  low  value. 

The  time  required  for  complete  vaporization  of 

magnesium  particles  of  size  400  pm  was  chosen  as 

0.1  seconds,  a compromise  between  single  particle 

• 2 1 

burning  experiments  which  indicate  longer  times 
and  observations  of  the  position  of  maximum  in- 
tensity14, which  indicate  shorter  times.  In  the 
case  of  finer  particles  a burning  time  proportional 
to  the  square  of  the  initial  diameter  is  to  be 
assumed,  in  accord  with  both  theory  and  experi- 
ment.21 

For  the  graybody  radiation  constant,  ji,  the 
value  2.6  x 10'11  cal/(sec  g deg4)  was  chosen  as 
giving  a reasonable  visible  plume  length  with 
the  plume  angle  chosen. 

There  is  some  experimental  evidence~~  for 
emissivities  as  high  as  0.5  for  MgO  and  Al^ 
particles  at  high  temperatures.  This  would  result 
in  nearly  blackbody  radiation  intensity  from  the 

envelope  of  the  smoke  cloud. 

From  the  value  of  a^  which  was  found  to  give  the 
best  results  in  the  present  calculations  the  gray- 
body radiation  efficiency  of  the  smoke  cloud 
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envelope  is  found  to  be  much  less  than  unity,  per- 
haps about  0.1.  This  implies  a much  lower  particle 
emissivity  than  in  the  experiments  referred  to.23 
In  support  of  the  present  results,  if  the  smoke 
were  relatively  black  one  might  expect  the  atomic 
and  molecular  spectral  features  to  be  relatively 
swamped  by  a large  graybody  continuum,  contra.ry  to 
what  is  observed.1'19 

The  sodium  D radiation  was  computed  theoreti- 
cally. Various  approximations  made  in  this  calcu- 
lation probably  affect  the  absolute  magnitude  of 
the  computed  radiation  more  than  they  affect  its 
dependence  on  temperature  and  sodium  atom  concentra- 
tion. Therefore  an  adjustment  factor  was  added 
for  computing  the  absolute  magnitude.  The  value 
of  this  factor  is  closely  dependent  on  the  value 
chosen  for  the  graybody  constant  through  the  re- 
quirement mentioned  earlier,  that  the  sodium  D 
radiation  from  a Mk  45  Aircraft  Parachute  Flare 
equal  about  8%  of  the  total  heat  of  combustion. 

It  is  found  that  the  value  2 satisfies  this  con- 
dition well  enough  for  large  (11  cm)  flares.  This 
is  also  reasonably  close  to  the  hoped-for  value  of 
unity  which  would  have  been  obtained  if  the  theo- 
retical derivations  could  have  been  done  exactly. 

The  model  predicts  greater  outputs  for  smaller 
diameter  flares  than  for  large  ones.  This  is  con- 
trary to  experience,  and  may  be  caused  by  any  of 
a number  of  the  assumptions  in  the  model,  such  as 
the  assumption  of  complete  combustion,  the  manner 
of  computing  admixed  air.  or  the  approximations  in 
the  radiative  transfer  model  of  the  emerging 
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radiation.  Here  it  has  been  dealt  with  by  using  a 
lower  radiation  constant  for  smaller  flares--as 
small  as  1.2  for  a 3 cm  diameter  flare.  Since 
there  is  no  theoretical  rationale  for  doing  this, 
the  model  as  it  now  stands  should  not  be  used  to 
predict  the  result  of  varying  flare  diameter. 

The  ratio  1 /I  was  more  or  less  arbitrarily 
1 2 

taken  as  10.  This  value  is  reasonable  in  view  of 
a computation  of  flame  temperature  based  on  the  as- 
sumption that  the  intensity  maximum  equaled  the 
equilibrium  intensity  at  this  wavelength.13  Since 
the  result  was  in  the  vicinity  of  the  adiabatic 

temperature,  1 /I  must  be  much  greater  than  unity. 

1 2 

Fortunately,  the  model  predictions  turn  out  not  to 
be  critically  sensitive  to'  this  parameter. 

An  attempt  was  made  to"  determine  this  ratio 
from  the  broadened  line  spectrum  in  a systematic 
fashion,  using  the  ratio  of  the  apparent  width  at 
half  height  of  the  absorption  dip  at  line  center 
to  the  apparent  width  at  half  height  of  the 
entire  line.  However,  it  was  realized  that  the 
continuum  also  had  to  be  included  and  that  this 
greatly  complicated  the  computation.  Therefore, 
this  approach  was  abandoned. 

One  of  the  most  difficult  parameters  to  assign 
is  the  amount  of  admixed  air.  The  computed  results 
are  fairly  sensitive  to  this  factor.  As  it  is  in- 
creased the  optimum  fuel /oxi di ze r ratio  increases. 
For  reasons  given  earlier  the  proportion  of  air 
was  assumed  to  increase  linearly  with  distance. 

The  proportionality  constant  was  taken  such  as  to 
give  a ratio  of  two  parts  air  to  one  part  original 
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composition  by  a distance  200  cm  out  in  the  plume. 

It  is  possible  that  a larger  value,  say  3,  in  com- 
bination with  changes  in  the  other  parameters  might 
give  a better  overall  fit  of  the  available  data. 

Partly  on  empirical  and  partly  on  theoretical 
grounds  it  was  decided  that  this  number  should  be 
proportional  to  the  2/3  power  of  pressure  and  in- 
versely proportional  to  flare  diameter  for  other 
than  10.8  cm  flares  burned  at  one  atmosphere. 

Computed  and  Experimental  Values  for  Specific  Cases 

The  temperature  and  flame  brightness  along  the 
length  of  the  plume  were  computed  for  a flare  of 
composition  58%  magnesium,  37.5%  sodium  nitrate, 
and  4.5%  epoxy  binder,  and  of  diameter  10.8  cm, 
and  the  results  are  presented  in  Fig.  2.  Since 
most  of  the  parameters  were  adjusted  by  comparison 
with  observations  on  this  flare,  the  fit  of  com- 
puted to  experimental  results  would  be  expected 
to  be  good. 

Briefly,  the  length  of  the  visible  part  of  the 
plume  (to  a decrease  of  brightness  by  a factor  of 
ca.  106)  is  seen  to  be  about  two  meters.  The  so- 
dium D output  integrated  over  the  entire  plume 
length  is  8.7%  of  the  total  combustion  energy  and 
the  brightest  spot  is  about  20  cm  from  the  surface, 
all  in  accord  with  observation. 

Trials  with  various  proportions  of  admixed  air, 
and  with  other  parameters  constant,  are  shown  in 
Fig.  3a.  The  maximum  in  computed  temperature,  as 
well  as  the  integrated  sodium  D output,  seem  to 
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FIGURE  2.  Computed  flame  temperature  and  flame  brightness  along  the 
plume  axis  of  a flare  10.8  cm  diameter,  58%  Mg,  37.5%  NaNO^,  4.5% 
epoxy  binder.  This  is  the  plot  from  the  run  listed  in  Table  A-3. 
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come  at  an  ai r/ composition  weight  ratio  (200  cm) 
between  3 and  4.  The  value  2 gives  nearly  the  max- 
imum output,  and  as  mentioned  earlier  other  fac- 
tors besides  maximum  output  were  considered  in  the 
selection  of  this  ratio. 


It  is  interesting  to  see  what  temperature 
could  be  expected  for  various  proportions  of  air  if 
no  heat  were  lost  from  the  system.  For  this  pur- 
pose we  have  computed  adirbatic  flame  temperatures 
as  a function  of  degree  of  vaporization  of  mag- 
nesium and  of  proportion  of  air,  for  three  differ- 
ent compositions  of  a Mg/NaNO  /binder  flare.24 
These  temperatures  are  presented  as  topographic 
surfaces  in  Figs.  3b,  c,  and  d.  Composition  paths 
representing  the  various  proportions  of  air  shown 
in  Fig.  3a  are  superimposed  on  the  appropriate 
adiabatic  surface,  Fig.  3c. 

A comparison  of  experi mental  and  calculated 
luminous  efficiencies  at  various  atmospheric  pres- 
sures and  for  a number  of  fuel /oxi  dizer/bi  nder 


ratios  is  presented  in  Table  1. 

The  agreement  between  prediction  and  experiment 
is  fairly  good,  considering  that  the  experimental 
values  came  from  three  different  groups  of  workers. 
Generally  it  can  be  said  that  the  calculation  pre- 
dicts a greater  effect  of  pressure  on  fuel-rich 
mixtures,  and  a lesser  effect  of  pressure  on 
stoichiometric  or  oxidizer-rich  mixtures  than  is 


found  experimentally. 

A comparison  of  experimental 26  and  calculated 
luminous  efficiencies  in  three  atmospheric  en- 
vi ronments--ai r , nitrogen,  and  argon--is  presented 
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TABLE  1.  EFFECT  OF  AhffilENT  PRESSURE.  Coroarison  of  predicted  and  observed  luminous  efficiencies,  cd 
sec/g,  for  magnesium/sodium  nitrate/binder  flares.  From  left  to  right  the  mixtures  are  Increasingly 
fuel -rich. 
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in  Table  2.  The  predicted  values  are  seen  to  fol- 
low the  experimental  values  semi quantitati vely  as 
fuel  percent,  binder  content  and  the  environment 
are  changed. 

The  biggest  disagreement  is  that  computed 
values  for  a fuel-rich  mixture  in  a nitrogen  en- 
vironment are  much  lower  than  the  experimental  ones. 
This  is  similar  to  the  discrepancy  for  fuel-rich 
flares  at  low  pressure  (Table  1).  Also,  the  com- 
puted superiority  of  argon  over  nitrogen  is  greater 
than  that  found  experimentally. 

The  computed  effect  of  going  from  5%  to  8%  bin- 
der at  the  highest  fuel  percent  is  the  reverse  of 
that  found  in  these  experiments;  however,  general 
experience  tends  to  confirm  the  computations  that 
illumination  efficiency  decreases  with  increasing 
binder  content  when  other  factors  are  held  constant. 

Experiments  and  calculations  to  compare  the 
effectiveness  of  various  sodium-containing  oxi- 
dizers are  presented  in  Tables  3 and  4. 

Table  3 compares  calculated  maximum  tempera- 
tures and  luminous  efficiencies  for  magnesium- 
fueled  formulas  with  a variety  of  sodium-containing 
oxidizers.  The  same  grain  diameter,  burning  rate, 
binder  content,  and  magnesium/oxidizer  equivalence 
ratio  are  assumed  for  all  cases. 

Of  the  five  other  oxidizers,  only  sodium  super- 
oxide is  predicted  to  give  greater  luminous  ef- 
ficiency than  sodium  nitrate.  The  higher  tempera- 
ture predicted  for  the  sodium  perchlorate  formulation 
is  more  than  compensated  for  by  the  tieing-up  of  the 
atomic  sodium  as  sodium  chloride,  which  doesn't 
radi ate  wel 1 . 
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TABLE  2.  COMPUTED  AND  EXPERIMENTAL2®  LUMINOUS  EFFICIENCIES  IN  VARIOUS  ENVIRONMENTS 
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TABLE  4.  EXPERIMENTAL  VERSUS  COMPUTED  LUMINOUS  EFFICIENCIES  FOR 
VARIOUS  FORMULAS  WITH  A SODIUM-CONTAINING  OXIDIZER 
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A comparison  of  predicted  luminous  efficien- 
cies with  experimental  results  for  various  sodium- 
containing  oxidizers  is  given  in  Table  4.  Each 
type  of  oxidizer  is  paired  with  the  corresponding 
"control"  formula  using  sodium  nitrate.  The  calcu- 
lations assume  the  same  percentage  compositions 
and  grain  diameters  as  had  been  used  in  the  ex- 
periments. However,  a common  burning  rate  of 
0.25  cm/sec  was  used  except  for  the  sodium  super- 
oxide calculations  where  the  difference  in  the  ex- 
perimental burning  rates  between  the  odium  super- 
oxide formula  and  the  sodium  nitrate  control  was 
large  enough  to  make  a large  difference  in  the 
cal  cul  atioris . 

The  model  is  fairly  good  in  predicting  the 
result  of  substituting  sodium  perchlorate  or  sodium 
iodate  for  sodium  nitrate. 

The  difference  between  Tables  3 and  4 regarding 
the  predictions  for  sodium  nitrite  are  due  to  a 
difference  in  fuel  /oxi  dizer  ratio  between  the  sodium 
nitrite  and  the  sodium  nitrate  formulas  in  the  case 
of  Table  4.  Nevertheless,  the  model  fails  to  pre- 
dict the  moderate  decrease  for  sodium  nitrite 
actually  observed.  This  may  be  the  result  of  the 
basic  model  parameters  selected,  which  tend  to 
predict  too  high  values  for  fuel-rich  formulas 
and  too  low  values  for  stoichiometric  formulas. 

The  large  decrease  in  output  caused  by  substi- 
tuting sodium  superoxide  for  sodium  nitrate  is 
probably  caused  by  the  very  high  rate  of  burning 
which  the  superoxide  flares  had,  and  which  cannot 
presently  be  predicted  by  the  model.  An  excessive 
rate  of  burning  probably  causes  numerous  inefficiencies. 
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One  of  these  is  lack  of  time  for  radiation  before 
the  plume  gets  mixed  with  cold  ambient  air.  This 
effect  can  befitted  into  the  model.  The  calcula- 
tion employing  the  experimental  burning  rate  pre- 
dicts nearly  half  of  observed  decrease.  The  other 
effects  may  include  incomplete  combustion  of  the 
magnesium  particles. 

There  is  apparently  also  a failure  in  predict- 
ing a decrease  in  output  with  addition  of  sodium 
azide.  The  seriousness  of  the  failure  is  uncertain 
since  the  experimental  comparison  was  only  a visual 
estimate. 
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DISCUSSION 

Possible  Improvements 


The  model  predictions  of  luminous  efficiency 
correlate  moderately  well  with  experimental  results 
with  Mg/NaNO  flares  of  a wide  variety  of  fuel 

3 

percent*  a wide  range  of  burning  conditions,  and 
performed  by  a number  of  investigators.  The  big- 
gest discrepancies  are  probably  due  to  the  large 
fluctuations  typical  of  flare  performance.  These 
are  a reminder  of  the  importance  of  many  subtle 
factors,  such  as  oxidizer  particles  sizes,  loading 
pressures,  cracks  in  the  grain,  etc.,  which  would 
be  difficult  to  model.  There  are  even:  serious 
problems  in  getting  agreement  for  a given  set  of 
flares  observed  with  different  experimental  setups. 
The  most  that  a simple  thermodynamic  model  can 
do  is  predict  the  average  of  a large  number  of 
experiments . 

Some  of  the  systematic  disagreements  with  ex- 
periment, such  as  the  differences  in  the  effect 
of  decreased  ambient  pressure  on  compositions  with 
different  fuel /oxi  di  zer  ratios  might  perhaps  be 
removed  by  reoptimizing  the  input  parameters-- 
plume  angle,  radiation  constants,  proportion  of 
air,  etc. 

While  doing  this,  a few  refinements  might  be 
useful.  For  instance  the  graybody  constant  might 
be  made  proportional  to  the  weight  percent  of  con- 
densed phase  present  in  the  reacting  mixture.  The 
rate  of  air  mixing  might  be  made  proportional  to 
the  percent  of  volatile  material  (binder  and 
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oxidizer)  in  the  original  composition,  or  even  to 
the  reciprocal  of  the  instantaneous  density,  on 
the  grounds  that  it  is  the  rush  of  gas  that  en- 
trains the  environment. 

There  are  several  major  refinements  which 
would  be  satisfying  from  a theoretical  viewpoint, 
and  which  might  also  bring  better  agreement  with 
experiment.  First,  the  inclusion  of  radiative 
feedback  would  enable  at  least  an  estimation  of 
burning  rates.  A consideration  of  feedback  is 
also  necessary  to  the  inclusion  of  a particle 
burning  model,  where  the  rate  of  vaporization  of 
the  metal  particles  would  be  expressed  as  a func- 
tion of  local  temperature  and  oxygen  content. 

A consideration  of  the  effect  of  local  oxy- 
gen content  on  particle  burning  would  probably 
improve  predictions  of  performance  of  fuel-rich 
compositions  in  oxygen-poor  environments. 

Applications  and  Limitations 


The  model  as  set  up  is  somewhat  open  ended. 
There  is  no  sharp  line  between  formulations  and 
conditions  which  can  be  computed  and  those  which 
cannot.  At  the  one  extreme  are  systems  which  can 
be  computed  quantitatively  and  reliably  simply  by 
furnishing  input  data  as  derived  in  this  report 
on  punch  cards  to  the  computer  program  as  presented 
in  the  Appendix.  At  the  other  extreme  are  systems 
for  which  the  computer  program  presented  here  would 
be  a useful  starting  point;  but  would  require  major 
re vi si ons  . 
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In  between  are  cases  where  after  minor  additions 
or  deletions  of  statements  to  the  program,  and/or 
the  examination  of  experimental  data  to  obtain 
certain  input  data,  quantitatively  or  qualita- 
tively reliable  results  could  be  computed. 

Thus,  simply  by  supplying  the  standard  NASA 
reactant  cards17  one  may  estimate  sodium  D radi- 
ation in  absolute  units  for  any  formulation  where 
MgO  (or  a condensed  phase  of  similar  emissivity) 
accounts  for  most  of  the  radiative  output.  <or 
instance,  formulations  of  the  type  Mg/binder/X, 
where  X is  a sodium-containing  oxidizer  may  be 
compared.  This  may  be  done  for  any  desired  ambient 
pressure  or  percent  composition. 

With  the  same  ease  one  may  compare  sodium  D 
radiation  among  formulations,  all  of  which  have  a 
common  graybody  emitter  of  properties  different 
from  MgO.  For  instance,  the  above  comparison  could 
be  made  using  aluminum  instead  o'f  magnesium  as 
fuel.  The  absolute  values  of  the  results  would 
not  be  reliable,  however,  unless  one  comparsion 
with  experiment  using  aluminum  fuel  were  made. 

Comparative  estimates  of  other  radiative  pro- 
cesses besides  sodium  D emission  may  also  be  made 
by  inserting  the  appropriate  emission  formula  at 
the  place  in  the  program  where  radiation  is  calcu- 
lated. For  the  case  of  atomic  line  emission  or 
emission  from  an  electronically  excited  state  of 
a polyatomic  species  of  relatively  narrow  band 
width  the  rate  of  emission  would  be  proportional  to 

[X]e"hv/kT, 
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where  [X]  is  the  concentration  of  the  emitting 
species  as  computed  by  the  thermodynamic  calcula- 
tion, and  hv  is  the  energy  of  the  upper  level.  For 
instance,  a comparison  of  the  emission  from  the  % 
state  of  BaCl  at  5200  ± 100  A for  the  formulas 
Mg/Ba(N0  ) /C  Cl  versus  Mg/Ba(C10  ) may  be  pos- 

3 2 6 6 4 2 

sible  in  this  way. 

The  expression  used  for  sodium  D radiation 
could  readily  be  adapted  to  other  lines  broadened 
by  passage  through  a region  of  high  optical  density. 

More  difficult  is  the  estimation  of  radiation 
in  absolute  units,  as  would  be  necessary  for  a 
comparison  of  two  different  lines  in  the  same 
formula.  This  requires  that  an  adjustment  factor, 
which  represents  the  radiative  lifetime  of  the 
excited  state  and  compensates  for  other  approxima- 
tions, be  added  to  the  expression  given  above. 

This  is  most  easily  done  by  calibration  against 
power  spectra  already  available.  For  instance, 
to  compare  Na-D  against  2ir  BaCl  radiation  in  a 
formula  including  Na,  Ba,  and  Cl,  it  is  first 
necessary  to  "calibrate"  the  computation  of  the 
BaCl  emission  (calibration  of  Na-D  has  already 
been  done)  by  means  of  an  available  power  spectrum 
for  a formula  which  also  contains  Ba  and  Cl  (but 
not  necessarily  Na). 

Of  still  greater  difficulty  is  the  computation 
of  emission  from  electronic  states  where  the  band 
spread  is  not  small  compared  to  the  center  fre- 
quency. Not  much  consideration  has  been  given  in 
this  project  to  a derivation  of  the  appropriate 
radiation  formulas. 
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Of  greatest  difficulty  would  be  the  computa- 
tion of  lines  which  overlap  in  a spectral  region 
of  high  optical  density.  Derivation  of  the  appro- 
priate radiation  formulas  and  their  insertion  into 
the  program  would  probably  require  a major  new 
effort. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  model  of  flare  plume  combustion  and  radi- 
ation as  presented  here  has  been  shown  to  give  pre- 
dictions of  luminous  output  which  are  in  qualitative 
or  semiquantitati ve  agreement  with  experiments  on 
magnesium  flares  of  a wide  range  of  compositions, 
and  burning  under  a wide  range  of  environmental 
compositions  and  pressures. 

It  should  be  useful  for  predicting  radiant  or 
luminous  output  with  new  formulas,  and  for  opti- 
mizing the  ratios  of  the  various  components  in 
the  proposed  formulas  because  it  makes  use  of 
all  of  the  thermodynamic  information  of  the  com- 
ponents in  a systematic  fashion.  The  method  used 
herefore  has  been  to  compute  single  thermodynamic 
parameters,  such  as  heat  of  combustion,  the  adia- 
batic flame  temperature,  and  the  emitter  concen- 
trations and  to  use  these  as  a qualitative  guide 
in  estimating  output. 

It  is  recommended  that  theoretical  and  ex- 
perimental input  parameters  be  obtained  from  the 
literature  such  that  the  model  can  be  used  for 
predictions  of  output  of  colored  flares  and  of 
infrared  decoy  flares. 

It  is  also  recommended  that  the  model  be  re- 
fined by  (1)  reoptimizing  the  arbitrary  input 
parameters  (2)  including  a heat  feedback  term  in 
the  basic  equation  and  (3)  including  the  dependence 
of  the  rate  of  fuel  combustion  on  temperature  and 
oxygen  concentration. 
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APPENDIX  A 

THE  COMPUTER  PROGRAM 

The  basic  computer  program  is  that  of  Gordon 

ft 

and  McBride  for  computation  of  chemical  equili- 
brium. The  unneeded- subroutines  ROCKET,  RKTOUT, 
FROZEN,  SHCK  and  DETON  were  deleted.  A new  sub- 
routine, INCRMT,  was  added  to  add  increments  of 
the  surrounding  atmosphere  and  of  magnesium  vapor 
to  the  reaction  mixture  corresponding  to  incre- 
ments of  time  and  distance  from  the  burning  sur- 
face, and  to  compute  radiation  losses. 

A few  statements  were  added  to  subroutine  SEARCH 
to  determine  the  species  number  of  the  particulate 
fuel  . 

Subroutine  THERMP  was  modified  to  allow  infor- 
mation on  the  assumed  conditions  of  combustion  and 
a specification  of  desired  output  to  be  read  in. 

An  entry  point  in  subroutine  EQLBRM  was  ar- 
ranged so  that  computations  would  loop  u'-'und 
through  INCRMT  and  EQLBRM  until  the  computed  flame 
temperature  had  dropped  to  a suitable  value  (now 
set  at  1000  K)  or  until  a limited  number  of  points 
had  been  computed. 

The  flow  charts  of  these  changes  are  shown  in 
Figs.  A-l  to  A-3.  A printout  of  the  modified  portions 
of  ti.e  NASA  program  is  given  in  Table  A-l.  A list 


S.  Gordon  and  B.  0.  McBride,  Computer  Program  for 
Calculation  of  Complex  Chemical  Equilibrium  Composi- 
tions j 'Rocket  Performance , Incident  and  Reflected 
Shooks j and  Chapman- Jouguet  Detonations , NASA  SP-273, 
Lewis  Research  Center  (1571).  N71-37775,  AD  727607. 
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TABLE  A-l  . 


IMPORTANT 

Name 

(dimensions) 

BD(5) 

BL(5 ,2) 

B LI (5 ) 

B L 2 ( 5 ) 

BRT 

BRIT 

B U F ( 6 3 ) 

CBB 

CNAD 

CON 

CON  C ( 10) 

DEN 

DTL 

DX 


VARIABLES  ADDED  TO  THE  NASA  PROGRAM 


Description 

Graybody  energy  radiated  since  the 
last  point  in  the  wavelength  band 
designated  by  the  corresponding  BL , 
calories/R/gram  of  initial  mixture. 

Wavelength  (ym)  limits  of  graybody 
radiation  band  desired. 

BL(  1 ,1 ) . 

BL ( 1 ,2 ) . 

Rate  of  recession  of  flare  surface, 
cm/sec . 

Flame  brightness  in  the  vicinity  of 
the  point  computed,  candela/cm2. 

Buffer  used  by  the  plot  subroutine. 

Adjustable  graybody  radiation  constant 
in  calories/R/deg4/sec/gram  of  ini- 
tial mixture. 

Adjustable  sodium-D  radiation  constant, 
di mensionless  . 

The  right  side  of  Eq . (1). 

g-moles  of  j-th  species/g  of  mixture 
including  admixed  environment,  and 
excluding  unvaporized  fuel. 

Density  of  flare  composition,  g/cm3. 

Sought-for  temperature  change  be- 
tween points. 

Current  distance  between  points. 
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TABLE  A-l  (cont.) 


DZRO 

EFF 

EMARAY(202) 

FLOWZ 

HSUBO 

KB 

KSP 

KRON 
LI  2 
NA 


Di  ameter  of  fl  are*  cm. 

Computed  graybody  emission  relative 
to  the  blackbody  radiation  from  a 
plume  of  the  same  diameter. 

Values  of  flame  brightness  stored 
for  plot. 

Flow  of  material  across  the  plane  of 
the  flare  surface,  g/sec. 

Enthalpy  of  reacting  mixture,  in-= 
eluding  environment,  excluding  unva- 
porized fuel  cal/R/g,  as  In  original 
NASA  program. 

Number  of  radiation  bands  desired. 

Number  of  species  concentrations  to 
be  printed  out. 

This  variable  equated  to  zero  to  com- 
pute adiabatic  temperatures.  Other- 
wise, equated  to  unity. 

1 1 /l s , the  effective  ratio  of  hot 
path  length  to  cold  path  length 
for  the  emerging  radiation. 

Species  number  of  atomic  sodium. 


NCF 


Condensed  species  number  of  the  fuel. 


NEND 


Maximum  number  of  points  desired. 


NMG  Species  number  of  the  fuel. 

NPL PTS ( 202 ) Number  of  points  to  be  plotted. 


NX 


Current  point  number. 


PLOTT 


Logical  variable.  "True"  if  a plot 
is  desired. 
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TABLE  A-l  (cont.) 


PLT 

Suppresses  plotting  of  the  first 
poi nt . 

POWR 

The  amount  of  admixed  air  is  propor- 
tional to  the  distance  to  a power. 
Has  generally  been  equated  to  unity. 

RADI 

The  integral  in  Eq.  (6). 

RBB ( 5 ) 

Graybody  radiation  for  current  point, 
cal/R/g  of  original  composition. 

RH  0(5) 

Mixture  density  including  unvaporized 
fuel , g/cm3 , 

RNAD(5) 

Sodium-D  radiation  for  current  point, 
cal/R/g  of  original  composition.. 

SBD(  5 ) 

Total  graybody  radiation  in  a given 
wavelength  band,  units  of  calories/g 
for  final  printout. 

SBL 

Total  luminous  output  due  to  gray- 
body radiation,  cd  sec/g. 

SPC(9 ,2) 

Names  of  species  whose  concentrations 
are  desired  in  the  printout. 

SRBB (5 ) 

Total  graybody  radiation,  units  of 
cal/g  in  final  printout. 

S RNAD ( 5 ) 

Total  scdium-D  radiatioi.  units  of 
cal/g  in  final  printout 

SRNADL 

Total  luminous  output  due  to  sodium-D 
line,  cd  sec/g. 

STLUM 

Total  luminous  output,  cd  sec/g. 

T ( 5 ) 

Current  temperature. 

T2ARAY ( 202 ) 

Stored  values  of  temperature  for 
plot. 

THETA 

Assumed  half-apex-angle  of  plume, 
radi ans  . 
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TABLE  A-l  (cont.) 


b 


8 


I 


r 


T1 

T I ME ( 5 ) 
HUM 

TWF 

WF 


Guess  of  initial  temperature,  K. 

Current  time  of  flight,  sec. 

Luminous  output  at  current  point, 
cd  sec/g. 

Assumed  time  for  complete  vaporiza- 
tion of  fuel , sec. 

Final  weight  percentage  of  vaporized 
fuel,  usually  equal  to  the  percen- 
tage specified  on  the  corresponding 
"reactants"  card. 


WFA  Weight  ratio  of  admixed  environment 

to  composition  by  200  cm  from  the 
surface,  x 100. 

WI  Initial  weight  percentage  of  vapo- 

ri zed  fuel  . 


WT (14) 

xc 

XCARAY ( 202 ) 


Weight  percentages  of  reactants. 

Current  distance  from  the  flare 
s urf ace . 

Stored  values  of  XC. 


I 
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FIGURE  A-2.  Modifications  to  Subroutine  THERMP 


Unmodified  THERMP 
"TT=T(IT)H 


Initializations 


Read  conditions  for 
next  computation 


Runs 

completed 

(NEND=0)? 


Write  final 
point  no. 


Write  conditions 
and  headings 


Continue  with  unmodified 
THERMP  for  a final 
printout,  "IF(TT.NE.O)11 

r s 

( Return  ) 


Read  Species 


More  input 
(ROUTE^COND 
or  BLANK) 


Bands 
to  be 

^computed?. 


Read  and 
write  band 
limits. 


Write  column 
headings.  Compute 
initial  plume 
conditions 


Subroutine  EQLBRM.  \ 
Compute  and  write 
entire  reaction  path,/ 


Results  ^ 
to  be  plotted 
(PL0TT=T)  . 

\ / 


of  the  more  important  variables  introduced,  and 
their  definitions,  is  given  in  Table  A-2.  Table 
A-3  shows  a listing  of  input  and  output  from  an 
actual  run. 

Following  the  control  cards,  the  usual  cards 
for  the  NASA  program  must  be  furnished,  beginning 
with  the  THERMO  data  cards  if  the  thermodynamic 
data  is  not  furnished  on  tape.  The  format  is  un- 
changed from  that  described  by  Gordon  and  McBride, 
and  the  details  may  be  found  there  or  may  be  dis- 
cerned from  the  example  given  below. 

Because  of  the  modifications  there  are  a few 
special  restrictions  on  the  input  furnished  here: 

1.  The  first  reactant  card  must  be  for  the 
fuel,  assumed  to  be  an  element  in  particle  form. 

2.  The  last  reactant  card  is  for  the  in- 
tended ambient  environment.  All  elements  listed 
on  this  card  must  be  listed  on  one  of  the  previous 
cards  (in  order  to  avoid  a singular  matrix  at  the 
first  point),  though  this  need  be  only  in  a minute 
quantity.  For  example,  if  the  environment  contains 
argon,  one  of  the  cards  before  the  last  must  list 

a non-zero  amount  of  argon. 

3.  All  reactants  must  be  listed  as  fuels,  11 F". 

4.  On  the  "namelists"  card  only  the  "HP" 
option  is  to  be  used.  I0NS=T  may  not  be  used. 

" FPCT=T , MIX= 100 . " are  to  be  included. 

Immediately  following  the  "namelists"  input 
card  is  a card  with  the  word  "CONDitions" , fol- 
lowed by  a card  listing  the  assumed  combustion 
conditions--flare  diameter  and  burning  rate,  etc. 

These  conditions  and  the  input  formats  are  listed 
in  statements  numbered  32  and  31,  resp.  in 
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TABLE  A-2 

MODIFICATIONS  TO  NASA  PROGRAM 


Additions  or  changes  are  indicated  by  brackets  in  A throus.’  C. 


A.  Subroutine  SEARCH 

COMMON/ 1 NOX/'  !DEBUG#CONVG#  TP#HP#  SP#  ISV#  NPP#  MGLES#NP#NT#NPT#NLM. 

1 #N$#KMAT#  IMAT#  I01#N0F#N0MIT# IP#NEWR#NSUB#NSUP#  ITM#CPCVFR#CPCVEQ 
2 # IQNS#NC#NSERT# JSQL# JLIQ#KASE#NREAC# IC# JS1# VOL#SHOCK/ IT#NFZ/CALCH 

3 # IQS AVE#  LSAVE# I SUP# I SUB/ JTNUM#NX#NEND#  TWF  #Wl#WF#WFAiE$TR#BSTR 

4 /RATE/  DX#  TI  #KR0N#  DXF#  FLOWZk XZRO# A3#  WT(1<0#  AZRO#  L12 
a_J3#NA#NMG#NCF#lLB  ( 15 ) #NOUT  ( 1Q)#KSP,»BL  (3/2)  #KB#ILA(  15/  23#  2 } 

c conoense'o  species 

c 

NC*  NC+1 
TEMP (NC# 1)“  Tl 
TE.MP(HC#2)«  T2 

— IF  (MTU).NE,NAMEU#1)>  GO  TO  12 

11  IF  (MT(2).NE.IZERO.AND,MT'(2).NE.  JZEROI  GO  TO  12 
NCF  - NC 

NMG  ■ NS 

12  CONTINUE 


B.  Subroutine  EQLBRM 

1 SINGC  • .FALSE. 

C IF(NX.EQ.O)  CALL  JNCRMT 
PIE  • 0, 


43  IF  t.NOT.CONVG)  GO  TO  62 
Si  -I  i *=  CnN 
““  CALL  INCRMT 

IF (NX, EM. MEND)  GO  TO  1000 
IF(«X.GT. 10. ANO.TT.LT. 1000)  GO  TO  1000 
_ GO  TO  1 
62  CALL  MATRIX 
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TABLE  A*  (Coni) 
C.  Subroutine  THERMP 


Reproduced  from 
best  available  copy. 


COMMON/  'ISC/EWN, SUMN/TT/ SO/ ATOM! 3/ 1<H ) /LLMT ( 15 ), 00(15 )/BOP( 15*2) 

1 /T’,*/TLf''r./T,tI0/THIGH/PP/CPSUM/QF/EQRAT/FPCT/R/RR/H$UB0/AC(2)/AM(2). 

2 /HPP(2)/RH(2>/  VMINU)/VPL5(Z)/WP(2)/DATA(22)/NAME(15/5) 

3 /A '!,>■.(  l5/5)/P6CWT(l5)/ENTH(l5)/FAZU5)/RTEMP(15)/FOX(15)/DENS(l5) 
r A >RM(jP/RM;(15)/TLN/CR/QXF(15)/5NNL/ENSAVE/ENLSAV  /TRACE/$IZE/BD(3) 

COMMON/ )>„'X/  IDEBUG/CGNVG/TP/HP/SP/ ISV/  NPP/  M0LE5/NP/NT/NPT/NLM 

1 /N5/KMAT/ I M AT  > IQI/NGF/NGMIT/ IP/NEWR/NSUB/NSUP/ ITM/CPCVFR/CPCVEQ 

2 /IONS/i'C/NSERT/  JSOL/ Jl I Q/KASE/NREAC/ IC/ JS1/ VOL/ SHOCK/ IT/NFZ/CALCH 

[3  /IQSAVE/LSAVE/ISUP/ISUB/ITNUM/NX/NEND/TWF  /WI/WF/WFA/CNAO/  CBB 
A /POUR/  DX/  T1  /K*ON/  DU/  FIOWZ/XZRO/ A3/  WT(1A)/  DZRO/  L12 
5 / 1A  / ‘IMS#  NCF,  L L 8 ( 15)/  NOUT(lO)/  KSP/  BL(3/2)/  KB/  LLAU5/25/2) 
C0MM0N/0UPT/FMT<30)/FP(A)/FT(A)/FH(A)/FS(A)/FM(A)/Fv(A)/FD(A) 

* 1 /FC(A)/FG(4)/FB/FMT13/Fl/F2/F3/FA/F5/FL(A)/FMTI9/FA1/FA2 

2 /FR1/FC1/FN<A)/FR(A>/FA'<4)/FI<4)/FMT9X/F0 

EQUIVALENCE  (K/ISV  )/(VL/P)/(UV/HP)/(TP/TV)/(SP/SV) 

OATA  FUU/AHU/ -C / 

EOATA  CGND/ AHCOND//BAN0/AHBAN0// SPEC /4HSP EC// BLANK/ 1H  / 

OATA  SPCN(2)/AHNA  //OMIT/AHOMIT/ 

00  902  IT-l/NT 
TT  « KIT) 

BEG  INNING  OF  MODIFICATIONS, INITIALIZATIONS, 

00  215  I-l/NLM 
M * 0 

DO  213  J-l/NREAC 
DC  212  R«l/5 

IF  (NAME(J/K).NE.LLMT(I))  GO  TO  212 
M ■ M+l 

LLAU/M/l)  • J 
LLAU/M/2)  ■ K 
212  CONTINUE 
213  CONTINUE 
LL8(I)  « M 
215  CONTINUE 
READ  (5/31) 

>0  KSP  ■ 0 
K(-  - C 
00  52  1*1/5 
BTU)  ■ 0 
00  52  J*l/2 
52  BL ( 1/ d ) - 0 
00  56  I »i/9 
MOUT(I)  » 0 
DO  58  J«l/2 
58  SPC(I/J)  « 0 

I 

, READ  NEXT  SET  OF  BURNING  CONDITIONS, 

I READ  (5/31)  OZRO/DEN  /TI/BRT/THETA/WI/WF/  TWF/P0WR/WFA/CNAD/L12 
1 /CBB/  KRON/OX/DTL/NEND  /PLOTT 
I IF  (NEND.EQ.O)  GO  TO  AO 


TABLE  Ar2  (Cont,) 

C.  Subroutine  THERMP  (Cont.) 

'/'■  I T c ( t>>  3 2 > 

.-2  F c-'-hT  (6781  OZRO  T!  BRT  THETA  Wl  WF  TWF  POWR  WFA  CNAD  L 

11?  CBG  KROM  OX  OXF  NEND  PLOTT  } 

WRITE  (6/31)  DZRO/DfcN  /TI/BRT/THETA/WI/WF/  TWF/ POWR, WFA/ CNAD/ U.2 
1 /C»8>  KRON/DX/DTL/NENO  /PLOTT 

31  FORMAT  (1H  /F5 , 2/ F5 .2/ F5. 0/ F5.2/ F5,3/ FA . 1/ F5 . 1/ FA . 2/ FA . 1/ F6. 0/ FA.  1 
1 /I3/66. 1/12/FA. 1. FA, O/IA/12) 

50  READ  (5/5),}  ROUTE 

51  FORMAT  (AA) 

IF  (ROUT6„EQ»CQNO;OR. ROUTE. EQ, BLANK)  00  TO  79 
IF  (ROUTE, EQ. BAND)  CO  TO  53 

C READ  SPECIES  FQR  PRINTOUT  OF  EQUILIBRIUM  CONCENTRATIONS 
READ  (5/33) ( (SPC (l/JJ/J-1/2)/ 1-1/9) 

33  FORMAT  (18AA) 

00  37  1-1/9 

IF  (SPC(I/D)  38/  33/  3A 
3A  DO  36  J-l/150 

IF  (SUB(J/l).EQ.SPC(I/l).AN0,SUB(J/2),EQ.SPC(I/2))  GO  TO  35 

36  CONTINUE 

SPC ( 1/2 ) « OMIT 
35  NPUT(I)  • J 

37  CONTINUE 
KSP  ■ I 
CO  TO  50 

38  KSP  *1-1 
GO  TO  50 

C READ  WAVELENGTH  BANDS  TO  BE  CALCULATED, 

53  READ  (5/5A)  ((BL(l/J)/J«l/2)/I-l/5) 

5A  FORMAT  U0F5.2) 

WRITE (6/57)  (<BL(I/J)/J-l/2)/I-l/3) 

57  FORMAT  (1AH0  BAND  LIMITS  ///3(F9.2/AH  TO  /F5.2)) 

DO  55  1-1/5 

IF(BL(I/1).EQ,0)  CO  TO  56 

55  CONTINUE 
KB  » I 
GO  TO  50 

56  K?  - 1-1 
GO  TO  5^ 

79  If  (KSP+KR.GT.  9)  KSP-  9-KB 

C TABLE  wgAriNGS, 

WRITE  (ft/fiO) 

WKI fE  (6/61)  ( (SPC ( I/J )/ J-l/2)/ 1-1/ KSP) 

BO  FORMAT  ( 1Q5HQ  fiy  0 1 S T T!MC(«)  TEMP  RAD 
1 EFF  MOLE  FRACTIONS/’  BANo"  ENERGIES 

R1  FORMAT  ( 55H  (CM)  (K)  CAL/R 

I /BX/  1FAA) 

C MORE  INITIALIZATIONS, 

NA  « 0 

DO  72  I-l/NS 

IF  (SUR(I/1),EQ,SPCN(2))  NA-I 
. 72  CONTINUE 


LUM  BRIGHT  SODIUM 
) 

CDS/G  CD/CM2  CAL/R 
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TABLE  A-2 

C.  Subroutine  THERMP  (Cont.) 

W * 0 
'ILiISS  « 

Ot  53  I»1,'ilESS 
<>?  w«w+P6C  ’T(  I ) 

DC  94  I *2/ ML ESS 
54  WT ( I ) * 100*PECWT( I )/••< 

XZRO  ■ DZFQ/ 2/TAN (THETA) 

FLOWZ  » .785398*02 RO**2*BRT*DEN 
NX  > 0 

CALL  EGL8RM 

IP( .NOT.PLOTT)  NPlPTSbO 
IP( .N3T.PL0TT)  GO  TO  30 
■ CALL  PLOTS ( RUF/63/  20) 

CALL  PACT0R(2.) 

EMARAY(l)  ■ EMARAYU) 

CALL  Pl0T(0./-U,;~3) 

CALL  PL0T(1.5#9./-3) 

CALL  SCALe(T2ARAY,5.,NPLPTS,l) 

CALL  AxjS(0,,0.,l9HTEMPERATURE,  BEGS  K,  19,5.,0.,T2ARAY<NPLPTS 
1 +l),T2ARAY(NPLPTS+2) ) 

. CALL  SCALE(XCARAY;7,,NPLPTS*1) 

1 CALL  AXISIO./O./  12HD I STANCE*  CM,-i2/7,,-90.,XCARAY(NPL'PTS+l  >, 

■ 1 XCARAYJNPLPTS+2) ) 

I CALL  SCALE(EMARAY,5.,NPLPT$,1> 

^ CALL  AXIS(0.,-7.,27H  LOGILUMINANCe)*  CO/SQ-CM  27,03. ,0./ 

, 1 EHARAY(NPLPTS+l)iEMARAY(NPLPTS+2)) 

1 XCARAY(NPLPTS*2)«-XCARAY(NPLPTS+2) 

| CALL  L1NE(T2ARAY,XCARAY,NPLPTS, 1/0,0) 

CALL  LINE(£MARAY,XCARAY,NPLPTS, 1,0,0) 

CALL  PL0T(7./-10.;-3) 

PL  T*0  # 

NPLPTSbO 
GC  TO  30 

40  'WRITE  (6/41)  NX 

41  FORMAT  (15) 

C CIV-TINjS  “ITH  UNMODIFIED  SUBROUTINE 
IF (TT.NE.C. ) GG  TO  000 
IF(NPT.fia,0)  GO  TO  1000 

e<  o k • o 
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TABLE  A-2  (Cont.) 
D.  Subroutine  INCRMT 


S' '•'•ROUT  INF  INCRMT 

Jl  FV>  JIA  t-(  5)/T(  5 )/TlME{  5 )/RNA0(  3 )/RI)B  ( 3 )/  C0NC(10) 

1 / l L 1 ( 5 ) / 8 L2  ( 5 > * RHO  ( 3 > / SRNAO  ( 5 ) , SRBB  ( 5 ) / GMR  ( 1 A ) 

2 / $:•  l‘.!(  B8E(l8/2)/  BU(ll) 

UCOlCftL  HP/SP/TP/  CONVG/NEWR/ IONS; MOLES; EQL/  FRQZ/ VOL 

Cr-  •'-..•/i'U’TC/'  PtPTSjXCARAYt 202 )/ T2ARAY (202)/EMARAY{ 202 } /8UF(fc3) 

CGMKU'l/PLTR/  PLT 

C''!*',ON/  ISC/ENN/$UMN/TT  / SO/ ATOM ( 3/  101 ) / LLMT (15)jB0(15)/B0P(15/2) 

1 /TM/TLQW/YMID/TH!GH/PP/CPSUM/QF/EQAAT/FPCT/R/RR/HSUB0/AC(2)/AM(2) 

2 /HPP(2)/RH(2)/  VMIN(2)/VPLS(2)/WP(2)/0ATA(22)/NAME(15/5) 

3 /ANUM(l5/5)/PECWT<l5)/ENTH(l5)/FAZ<l3)/RTEMP<15)/FQX'(15)/DEN$<l'': 
A ,RH0P,RMW(15)/TLN/CR/0XF<15)/ENNL/ENSAVE/ENLSAV  jTRACE/S1ZE/BD(5) 

COMMON/ INOX/  IDEBUG/CONVG, fP/HP/SP/ISV/  NPP,  MOLES/ NP/ NT/ NPT/NLM 

1 /NS/KMAT/ IMAT/ IQl/NOF/NOMlT/ IP/NEWR/NSUB/NSUP/ ITM/CPCVFR/CPCVEQ 

2 ,I0NS,NC/NSERT/JSQL/JLIQ/KASE/NREAC/IC/JS1/ VOL/SHOCK/ IT/NFZ/CALCH 

3 > IQSAVEj LSAVE; I5UP/ ISUB/ ITNUM..NX/NENO/TWF  ,WI/WF/WPAiCNAD/  CBB 
. A /POWR/  OX/  TI  /KRON/  DTI/  FLOWZ/XZRO/ A3/  WT(14)/  DZRO/  L12 

5 , NA / NMG/  NCF/  LLB(15)/  NOUT(lO),  KSP/  BL<9>2)>  KB/  LLAJ15/25/2) 
COMMON/ SP EC ES /COE F (2/7/15Q)jS(150)/EN(150/13)/ENLN( 150) /HO ( ISO ) 
l/DELN<L50)/A(15/150)/SUB(150/3)/lUSE(150)/TEMP(50/2)/SLN(130j 
DATA  BBE  / 1800,/ 1900,/ 2000,/ 2100,/ £200 • 

1 / 2300, , 2400. /2500./2&00,/ 2700,/ 2A00,/ 2900,/ 3000,/ 3100./ 1200, 

2 / 3 300,/ 3400./ 3500,/  ,001/ ,0017/ ,00255/  .0034,  ,0052/  ,007/  ,009 

3 / .012/  .015/  ,018/  ,0215/  .0253/  .03/  .0345/  .0395/  .044/  .049 
h / .035  / 

OP TF(A/8)«(( A* 1, )**{"!,/{ A*B*B )))*(li*((A+l»)**(-l,/(B*B)))) 

T(  1 ) * TT 

IF (NX, E 0.0)  TT»TI 

C RELOCATION  OF  STORED  QUANTITIES 
30  2 1*1/4 

T (6-1)  • T (5-1) 

TIME (6- 1 ) « TIME ( 5-1 ) 

RNA0(6-I)  « RNAO ( 5-1 ) 

RBB  (6-1)  * RBB  (5-1) 

SRNAO (6- 1 ) * SRNA0(5-I) 

3188(6-1)  » SP.BB  ( 5-1 ) 

HMC  (6-1)  * RPQ  (5-1) 

01  AM (6-1 ) « 01 AM ( 5-1 ) 

2 CnNT  I r 

C Flame  ?" I iHTNESS/  SYSTEM  DENSITY/  MOMENTUM  CONSERVATION/  GRAY 
C Q00Y  EFFICIENCY. 

IP  (NX-1)  32/  9/  7 

7 BP  I T ■ 10O,*TLUM*HHO(3>*DIAM(2)/(100+WT(NREAC))/(TIME'(2>-TIME(3)> 

P RH.M2)  • ,01219*PP/ENN/TT*( 100+WT (NREAC ) > /W 
IP  (NX.LT.2)  00  Tn  h2 
Rh02  « PP*RMW  ( NRE  AC  ) / R/T,(  2 ) 

Ci-N  • ( RH02+RH0  ( 2 ) ) /RHU*WPA*POWR*  (XC/200 ) *+ ( POWR-1 ) /2E4- ( 1*WT 
1 (NKkAC)/100)*2/(XC+XZKU) 

IF  (CO’j.GT  ,0)  WRITE  (6/84) 

04  FCRI'.AT  (3AH  FORWARD  MOMENTUM  CHANGE  POSITIVE  ) 

EFF  » 3,64E13*RBB(2)*RH0(2)*0IAM(2)/TT**4/(TIME(2)-TIME(3)>/ 

1 ( 100+wT (NREAC  ) ) 

C PRINTOUT  OF  RESULTS  FOR  LAST  POINT 
00  9 I-l/KSP  *■,. 


TABLE  A-2  (Con*.) 
D.  Subroutine  INCRMT  (Cont.) 


•ii  = \CjTU) 

9 C <•(!)  = l) 

??  ,vMT6  (6#*3)  <“X#XC..  T IMb  (2 ).»  T(2)#  RBB(2)#  TLUM,  eRlT#  RNAD(2)#EFF 
1 , ( CC.  C ( I ) # I «1#RSP  )#  <HD(I)#I*1#K6> 

53  F‘  *'(AT  QMO/IS#  F6.1#  F6.4#  F7.1#  F6.1#  FB.l#  E9.2#  F7.3#  F5.3 
1 #9F8. 5 ) 

IF(PLT  ,EQ.O, ) GO  TCI  42 
I!PL«>TS«'IPLPTS  + 1 
XCARAY(NPLPTS)»XC 
T2ARAY ( MPLPTS  J *T( 2 ) 

E"ARAY(-»pLpTS)  « ALOGiO(BRIT) 

IF ( 9P.IT. LE . .0001  j EKARAY(NPLPTS)«-4, 

*2  PLT«1. 

c KF.ADJUSTMI-NT  OF  INTERVAL  between  points. 

IF  (NX.LTj.6)  GO  TO  12 

ROT  • A3S(T(5)/6+T(4)/3+T(3)/2-T(2))|i,.6  /DTL 
OX  ■ DX*(l.3+,7*R0T)/(I+R0T) 

12  NX  » Nx+1 

13  XC  » XC+DX 

:c  (Nx-i)  io#io#u 
\o  xc«o 
stluh  ■ 0 

SBL  « 0 
30  34  I«l#KB 
3 <•  SBO(I)  « 0 

c EXIT  on  completion  OF  calculations. 

11  IF(NX,EQ,NEND)  go  TO  221 

IFfNX.GT. 10. ANO.TT.LT, 1000)  GO  TO  221 

C PLUME  DIAMETER#  FLIGHT  TIME,  AND  WEIGHTS  OF  FUEL  AND  ENVIRONMENT 
C FOR  NEXT  POINT, 

OIA". ( 1 ) * DZ.RU*(Xr.  + XZPD)/XZRQ 
IF  (NX-1)  221#  15 i 16 
15  T I m g ( 1 ) « 0 
GO  70  17 

It  TIME(l)  ■ T IMF ( 2 ) +OX*RHO (2)*7B,54/FL0WZ/( 100+WT ( NREAC ) ) *DI AM( 1 ) **2 
17  WT ( 1 ) « WI+(WF-WI)*(TIME(1)/TWF) 

I F ( T(l).GT.WF)  v>T  ( 1 ) F 

wT ( i;Rt AC ) ■ WFA*( XC/200, )**PQWR 

»>  a C 

CC  20  I S l/NREAC 
20  Y>  * *T(  I )+W 

c heat  ccmtent  of  unvaporizeo  fuel 

HMu  • 0 
TF 5 ■ TT 
LM  t ' NMG 

IF  (TT.LT,TEMP(NCF#1))  LMG«LMG-1 
IF  (TT.C.T.TEMP(NCF#2) ) TEPkTEMP(NCF>2> 

00  30  I»i#5 
II  ■ 6~I 

30  H«G  • (HMG+C0EF(2#II#LMG)/II)*TEP 
HMG  » HMG+CQEF ( 2/6#  LMG)  A-15 


TABLE  A£  (Cont.) 

D.  Subroutine  INCRMKCont.) 

C I I T I -L  rUEGPATlJ..  OF  LINE  SHAPE, 

I'(  'X,*r.l)  GP  TO  3*> 

XL2=rL.=AT(Ll2> 

R.MH  a ’PTF < XL2>  *>.2) 

, X AV  = -.2 

n-  31  1=1^13 
A Av  = X.-JAV+.A 

31  RADI  = 0 AD1+4*QPTF ( XL2s  XWAV ) 

; X A Vs" 

DC  32  1«1>12 

I XaAV  » XwAV+,4 

32  RA01  ■ RA01+2*0PTF ( XL2/ XWAV) 

' RADI  « RADl*.2/3+AL0Gtll2+l.)/27.04 

RAOl  ■ 2*RA01 

SRNAO(l)  • 0 
SRBB(l)  ■ 0 
GO  TO  37 

C RADIATION  COMPUTATIONS  FOR  NEXT  POINT. 

35  IF  (NA  , EQ.O)  GO  TO  36 

RNAO(l)  • 3.75E5*W*EXP(-24416/TT>*{EN<NA,l>*ENN*TT**,5/DIAM(l> 

1 /AL0G(L12+1.)  )**'tg*RADl*CNAD*(TlME{l>-TlME(2)) 

SRNADU)  « SRNAD(2)+RNAOtl) 

36  ReB Cl ) • (TIME!  I ) -TIME  C 2)  )*CBB  *TT**4 
SRBB(l)  ■ SRBB ( 2 ) *RBB { 1 ) 

IF  (TT.GE.BBEI 1 > 1 j ) CO  TO  60 
BBL  » 0 
Gr’  TL‘  64 

60  DO  61  I«2il7 

IF  (TT.l6.88Et J/i>)  GO  TO  63 

61  CONTINUE 
WRITE  (6,62) 

62  FORMAT  (60H  TEMPERATURE  6XCEE0S  RANGE  OF  LUMINOUS  EFFICIENCY  TABLE 
1 ) 

63  P ■ { TT-BBE ( I j 1 ) ) / 100 

PA  ■ P+1 

Pl<  « P«1 

BBL  ■ (P*PM*88F(I-l/2)/2-PA*PM*BBE(I/2>+P|FPA+BBE(r+1^2)/2>*R8B 

1 (1) 

64  S8L  » SBL  +BBL 

IF  (KB, 60.0)  GO  TO  49 
OO  45  I»1jKB 

BU(ll)  ■ l,43BS325E4/8l(I*2)/TT 
BU(  1)  = 1.43B3325E4/8L(i>l)/TT 
DC  46  J»2,10 

46  8U(J)  » 8U( 1 )+• 1*1 J-l )*( BU( 11 )-BU( 1 ) ) 

31  * r>P(  1 ) **3  / { EXP  < BU(  1 ) )•»!  ) ■»'0U(ii)**3/CEXP(5uCIl ) )*•!) 

Of  47  j*2/10,2 

47  81  « B1+4*BU(J)**3/(EXP(9U(J) )-l) 

DP  48  J*3j  9t 2 

4P  81  » Bl+2*BU(J)**3/(EXP(BU(J>)-1) 

n 0 ( I ) » 91/3.*.15399MBUC2)-BU(l)  WT**4*CBB  *(T  IME  ( 2 )-TlME  ( D) 

+ 5 S P 0 ( I ) ■ SBO( I >+B0( I ) 

49  CONTINUE 

TLUM  = 33P*RNAD(1)+450*BBC 

STLUM  « STLUM+TLUM  A-16 
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Reproduced  from 
best  available  copy. 


C e\TnAL»Y  FUR  NEXT  POINT. 

37  riSUBO  ■ -.041U**NMG*(100+WT<NREAC)-W)-100.*{SRNADU)+SRBB(in* 

1 

0"  33  IsZiNREAf. 

38  HS.JBO  » HSU30<-WT(1)*ENTH(1)/RMW(I}/R 
hsubo  ■ Hsuao/w 

C OTHER  QUANTITIES  FOR  NEXT  POINT. 

HPP ( 2 ) «0, 

VPLS(2)»0. 

V«IN<  2 )»0. 

AM(2>«05 

00  220  \‘«1,NREAC 

?EC«T(n)  » WT(N)/W 

GMR ( N ) « PECWT(N)/RMW(N) 

HPP(2)»HPP(2)*  ENTH(N)*PECWT(N)/(RMW(N)*R) 

Am ( 2 ) • Am { 2 ) +PECWT ( N ) /RMW ( N ) 

220  CONTINUE 

AY  ( 2 ) ■ 1./AMI2) 

00  215  I«I,NLH 
80PU/2)  ■ 0 
IL  » UB(i) 

00  213  M«1,IL 
LLl  • LLA(I>M,1) 

L12  • LlA( I/M/2) 

80P ( l>  2)  « BOP ( I j 2 ) +GMR (LLl )*ANUM( LL 1#  LL2 ) 

213  CONTINUE 

30(1)  • BOP ( 1*2 ) 

215  CONTINUE 

NEWR  « .FALSE, 

93  RETURN 

C FINAL  PR  I* TOUT. 

221  SRNADL  « 330*SRNA0 ( l ) 

SPL  ■ A50*SBL 
SSNAu(l)  ■ RHSRMAO(l) 

SR3? ( 1 ) ■ R*SPBB(1) 

00  9u  I»liKii 

31.1(5)  * Rl  ( I / 1 ) 

■U  2 ( I ) = BL( 1,2) 

<50  seod)  = R*SBO { I ) 

ITS  (N,94>  SRNAO ( 1 ) > SRNAOL * SBL>  STiUM,  SRB6 ( 1 ) ^ ( BLl ( I ) j BL2 < I ) 
1 j SBO ( I ) > 1 s 1/KB  ) 

9A  F'-Rm,»T  (17H  SODIU'-O  ENERGY  F8.3>6H  CAL/G/21H0SO0IUM-0  LUMINOSITY 

1 F10.2j9h  CO-SEC/G  /21H  GRAY  BODY  LUMINOSITY  F10.2/9H  CD-SEC/G  / 

2 21H  TOTAL  LUMINOSITY  F10.2/9H  CO-SEC/G  /2AH0T0TAL  GRAY  BODY  E 
3NEKGY  F8,3>6H  CAL/G/2AH0BAND  GRAY  BODY  ENERGIES  /5(FA,1>1H-  FA.Ij 
A F10.3/6H  CAL/G/)) 

GO  TO  93 
EMO 
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TABLE  A-3.  INPUT  AND  OUTPUT  OF  A SAMPLE  RUN 


HACIANM 

Ml,» 

94. 

a.o 

f 1.74 

Nil.  N 1. 

O 3, 

77.9 

-111940.  S?M.19 

f ?.?4 

C 9. *5  H 4. 

14  Ol. 

.IS  N .1 

4.9 

•»M».  l?44.19 

r 

H I.H  0 *, 

0. 

a. a UMi.is 

r 

|NSf»t 

MOO(S) 

NATO  IS) 

0*IT 

CM 

MM 

CM? 

OMIT 

CMN 

C?0 

MM? 

CHS 

OMIT 

C?M 

c 

CM 

c? 

OMIT 

CH?0 

CM? 

C?H? 

C?M4 

OMIT 

ct* 

C3 

C* 

C9 

M4MCUSTS 

SlMMT?  MMuT 

. M 1. 

• raetiTi  mu  • loo. 

t 

CONOlTlOMi 

la.*  1.7  iota,  »?9 

•3  ?.  91.  .1  1. 

?tt.  ?• 

la  i.tc-u  i .? 

ii.  laa  t 

*MCICS 

CIS)  CO 

CO? 

M?0  M4  MOO 

•A*>S. 

US  ?»9  2.5 

3.9 

MAC 

M4CTANU 
HO  1.0000 

-0,0000 

-o.ocoo 

-0,0000 

-0,0000 

41,000000 

0.00 

s 

241.150 

f 

1.74000 

NA  1.0CC0 

l.COOO  0 

1.0900 

.0,0000 

-0.0000 

>7,900000 

-111140,00 

s 

244,110 

A 

2.24000 

c s.aioo  m 

I.H'Q  0 

1.113C  M 

• sooo 

-0,0000 

4.100000 

-11400,00 

l 

244.110 

a 

•0,00000 

h f.iaoo  c 

it*  ot'0 

-o.ccoo 

-0.0000 

-0,0000 

0.000000 

0,00 

0 

244,110 

a 

•0.00000 

iNSitt 

M00IS) 

M420IS) 

OMIT 

c*« 

,.N 

CM? 

OMIT 

C 20 

Mm? 

Cm) 

OMIT 

CJH 

c 

CM 

c? 

OMIT 

tMja 

CN2 

C2M? 

C2M4 

OMIT 

C2N 

C) 

C4 

Cl 

NiMtitm 

140  CA40I  *r»TI  »MT|  A*  1.  , AACU1,  «|X  • 100.  4 

TOO  CCV  STARTS  UNSU»$C«1»T(0  A4AAV 

IAA0*  0H1  ClTfCUD  »Y  NAf|N»  AT  400MSS  COCllO 

CAlUO  4(Um  TA>«  AT  tlM  0177 

140  CA'Cl  4UAT2  *A*T,  »*•  i AACUT,  -IX  • 100.  > 

too  Mw  constants  ac*  u»i$uuc*  .»»io 

I4»0»  NuHif*  0044  CITICTI0  I f K*m|||*  4T  400*115  000110 
CAIK0  MOM  UN*  4T  UM  1*177 


•JK4TJ 


KASt 

0* 

T 

0.0,  C,0«  0.0,  0.0,  0,0.  0f0t  0,0.  0,0,  0,0.  0,0a  0,0# 

0.0»  0.0*  0.0,  C.O,  0.0,  0.0,  0.0,  0.0,  0.0,  0,0, 

0.0, 

0,0, 

0.0, 

0.0, 

o.o. 

a 

0.IU51,  0,0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0,0,  0,0,  0.0, 

O.C,  0.0,  0.0,  0,0,  0.0,  0‘,0,  0,0,  0.0,  0.0,  0,0, 

0.0, 

0.0, 

0.0, 

0.0, 

0.0, 

MU 

a. 

WMMC 

a, 

NSOM 

a. 

V 

0.0,  0.0,  0.0,  0.9,  0.0,  0*.0,  0.0,  0.0,  0,0,  O.O.o  0.0, 

0.0, 

0.0, 

AmO 

0. 11*01,  0,0,  0.0,  0.0,  0.0,*  0,0,  0.0,  0.0,  0,0,  0,0,  0.0, 
0.0,  0.0,  0.0,  0,9,  0,0,  0,0,  0,0,  0.0,  0,0*  0,0, 

0.0, 

0.0, 

0,0, 

0.0, 

0.0, 

(MTIO 

a. 

01 

a. 

aact 

t. 

n 

a. 

MU 

0,11*01,  0.0,  0.0,  0,0,  0.0,  0,0,  0.0,  0,0,  0.0,  0,0,  0.0, 

0.0, 

0,0, 

0.0, 

0.0, 

TA 

a. 

MA 

T, 

la 

a, 

TV 

a, 

UV 

a. 

IV 

a, 

ART 

a. 

smock 

a, 

oitn 

a. 

A- 1 8 


r ' ; ^>  ivr*.i  ' ^v>  v ’ » io* '*•' 


TABLE  A-3  (cont.) 


21 

4.1 

.0140 

14*1.4 

.4 

1.4 

•Ilf *02 

•0(4  .044  O.Q0000 

•0000) 

.00*2* 

.00042 

.0929) 

•09000  .199)2  .07*44 

22 

4.4 

.0202 

2041.4 

.3 

2.4 

.141*02 

■0(4  .042  0. 00000 

•0000* 

.00*1* 

•00074 

•002** 

•9000  0**32  «0D2 

21 

7.2 

•Ml* 

20«4.3 

.3 

1.4 

.241*02 

.C«9  .041  p.OOOCO 

.00007 

•00*0* 

•00094 

.00297 

•00090  ,19)19  ,09*77 

24 

7.7 

.0224 

2134.7 

.4 

3.) 

.341*02 

,01)  .040  7*00090 

•0001 t 

•00)9) 

•OOlll 

•00)27 

•9600  00309  *0299 

25 

4.2 

.0214 

22-7.4 

.* 

7.* 

.4(1*02 

.919  .09  1.00000 

• 0000 

.00)11 

.00123 

•00)3* 

•0090  014)1  .090)4 

24 

4.7 

.0247  2234*2 

.7 

10.4 

.*91*02 

.027  .0)1  4.00000 

•0002) 

•com 

•001)0 

•00)92 

•000  04772  0OTC7 

27 

•.2 

.3233 

2)04.4 

,1 

14.3 

•*9f *02 

.0)7  .0)7  9.00000 

•000)2 

•C0131 

•00149 

•00409 

•OOOOO  07)92  .1150 

2i 

4.* 

.0271 

2)54.2 

.9 

19.3 

.111*0) 

.030  .034  ft. 90000 

.0004) 

.001)4 

.00137 

•00429 

•00009  00201  .12912 

24 

10.1 

.0241 
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TABLE  A- 3 (cont.) 
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subroutine  THERMP,  and  may  also  be  seen  in  the  sam- 
ple output. 

At  this  point  a card  with  the  word  "SPECies" 
may  be  placed,  followed  by  a card  listing  up  to  9 
species  (less  the  number  of  desired  bands)  for 
which  mole  fractions  are  to  be  printed  at  each 
point  computed. 

Before  or  after  the  species  cards  one  may  place 
a card  with  the  word  "BANDs",  followed  by  a card 
listing  the  wavelength  bands  (in  tun)  in  which  it 
is  desired  to  know  the  amount  of  graybody  radiation. 
Up  to  five  bands  may  be  specified.  Formats  for 
species  and  bands  cards  are  given  in  statements 
numbered  33  and  54,  respectively  in  subroutine 
THERMP. 

Upon  encountering  a blank  card  or  a second 
"CONDitions"  card  the  program  proceeds  with  actual 
computations.  A second  blank  card  here  causes  the 
program  to  exit  after  the  computations  are  complete. 
Otherwise  more  input  parameters  may  be  furnished 
for  another  calculation  using  the  .same  fonmul  a 
and  percentage  amounts,  but  different  burning 
condi ti ons . 


